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I.  INTRODUCTION 


The  primary  interest  of  any  aviator  engaged  in  flight  is  control  of 
his  aircraft.  Without  this,  it  becomes  impossible  tc  maintain  the  air¬ 
craft  in  a  desired  flight  regime  and  consequently  the  flight  is  rather 
abruptly  terminated.  With  adequate  control  over  his  vehicle  however, 
the  aviator  is  capable  of  successfully  performing  many  various  tasks  to 
include  a  safe  final  landing. 

Historically,  the  development  of  the  helicopter  was  hampered  during 
its  early  years  of  evolution  because  of  the  lack  of  understanding  of  the 
factors  which  caused  stability  (or  instability)  and  therefore  adequate 
control  over  the  vehicle  was  difficult  to  achieve.  This  situation  was 
exactly  similar  to  that  of  the  early  development  of  the  airplane.  It 
was  only  when  the  principles  of  stability  and  control  could  be  under¬ 
stood  that  aircraft  could  be  developed  which  would  fly  as  the  designer 
intended  them  to  fly  and  that  helicopter  development  could  progress. 

Before  the  details  of  controlling  the  aircraft  can  be  fully  worked 
out,  some  understanding  of  the  aircraft's  inherent  stability  must  be 
attained  because  these  two  factors,  stability  and  control,  are  closely 
related.  Together  they  determine  the  flying  qualities,  or  handling 
qualities  of  the  aircraft.  The  aircraft  must  have  sufficient  stability 
to  maintain  a  certain  desired  condition  of  flight  and  to  recover  normal¬ 
ly  from  disturbing  influences  (wind  gusts,  for  instance).  Pilot  work¬ 
load  is  also  a  function  of  stability.  Since  adequate  maneuverabi  lity  is 
a  necessity,  the  aircraft  must  also  respond  properly  to  the  pilot's 
i nputs. 
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An  expansion  of  the  concepts  behind  aircraft  stability  will  be  dis¬ 
cussed  here  so  that  they  can  be  more  fully  understood  and  easily  dealt 


with  when  moving  into  the  area  of  aircraft  control. 

A.  EQUILIBRIUM 

A  helicopter  is  in  a  state  of  equilibrium  when  the  vector  sum  of  all 
forces  and  moments  on  it  are  equal  to  zero.  While  in  equilibrium  the 
aircraft  will  not  have  any  tendencies  to  accelerate  in  either  transla¬ 
tional  (no  unbalanced  forces)  or  rotational  (no  unbalanced  moments) 
directions.  Thus  the  aircraft  will  remain  in  a  steady  flight  condition. 

If,  however,  forces  or  moments  are  introduced  to  upset  this  balanced 
condition  (via  cockpit  control  inputs  or  wind  gusts,  for  example)  the 
helicopter  will  experience  an  acceleration  in  the  direction  of  the  un¬ 
balanced  forces  and/or  moments.  As  can  be  expected  from  Newton's  Second 
Law,  linear  accelerations  are  proportional  to  the  magnitude  of  the  un¬ 
balanced  forces  and  the  angular  accelerations  are  proportional  to  the 
unbalanced  moments. 

B.  STATIC  STABILITY 

The  static  stability  of  a  system  is  defined  by  the  initial  tendency 
of  the  system  to  return  to  equilibrium  conditions  following  some  distur¬ 
bance  from  equilibrium.  If  an  object  which  is  disturbed  from  equilib¬ 
rium  tends  to  return  to  equilibrium,  the  object  has  positive  static 
stability.  On  the  other  hand,  if  the  object,  upon  being  disturbed,  has 
a  tendency  to  continue  in  the  direction  of  the  disturbance,  then  the 
object  is  exhibiting  behavior  of  negative  static  stability.  Neutral 
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static  stability  exists  when  the  object  has  no  tendency  to  return  to 
equilibrium  or  to  continue  in  the  direction  of  the  disturbance. 


The 

Figure 


classic  physical  examples  of  static  stability  are 
-1  [Ref.  1].  Various  tendencies  of  motion  of  a  ball 


POSITIVE  STATIC  STABILITY 


EQUILIBRIUM  ENCOUNTERED 
AT  ANT  POINT  OP  DISPLACEMENT 


NEuTHAl  S'ATic  STABILITY 


Figure  1-1.  Static  Stability 


shown  in 
displaced 
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from  equilibrium  in  a  depression,  on  a  hill,  and  on  a  level  surface  are 
shown.  It  should  be  noted  that  positive  static  stability  is  the  desired 
response  in  most  situations. 

It  should  also  be  noted  that  there  are  quantitative  degrees  of 
static  stability.  This  has  to  do  with  the  forces  acting  on  an  object 
after  it  has  been  disturbed  from  the  equilibrium  condition.  An  example 
will  illustrate  this  point.  A  large  ball  in  a  shadow  depression  may 
have  a  force  of  one-half  pound  returning  it  to  the  equilibrium  position 
while  the  same  size  ball  in  a  very  steep-sided  depression  may  have  a 
restoring  force  of  ten  pounds.  While  both  these  systems  exhibit  posi¬ 
tive  static  stability,  the  second  is  more  positive  and  thus  a  more 
stable  system. 

If  a  stability  control  system  is  to  be  incorporated  into  the  heli¬ 
copter  then  the  magnitude  of  the  aircraft's  static  stability  terms  will 
be  one  clue  to  the  amount  of  force  the  control  system  must  have  to  be 
effective.  For  example,  if  a  system  has  negative  static  stability  in 
yaw,  some  control  feature  must  be  incorporated  to  allow  the  pilot  to 
keep  the  yawing  motion  under  control.  If  the  negative  yaw  static  sta¬ 
bility  term  is  of  large  magnitude,  then  the  control  system  will  have  to 
have  a  great  amount  of  power  to  control  the  tendency  of  the  aircraft  to 
diverge  in  yaw.  On  the  other  hand,  if  the  stability  term  is  only 
slightly  negative,  then  the  yawing  motion  will  not  diverge  as  quickly  as 
in  the  former  case  nor  will  the  controlling  power  needed  be  as  great  as 
in  the  former  case. 

For  a  complete  discussion  of  helicopter  control,  dynamic  stability 
must  also  be  included. 
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C.  DYNAMIC  STABILITY 


Dynamic  stability  refers  to  a  body's  resulting  motions  with  respect 
to  time  after  being  disturbed  from  equilibrium.  A  plot  of  displacement 
versus  time  will  reveal  the  dynamic  stability  tendencies  of  a  body.  All 
possible  responses  of  a  disturbed  body  can  be  seen  in  Figure  1-2.  Two 
general  modes  of  motion  exist,  oscillatory  and  non-oscil latory  (also 
called  periodic  and  aperiodic).  As  is  implied  by  the  term,  oscillatory, 
the  position  of  the  body  will  cycle  in  some  manner  about  the  equilibrium 
position. 

The  motion  of  both  the  oscillatory  and  non-oscillator  .odes  will 
also  depend  on  whether  or  not  such  motion  is  damped.  If  a  eventu¬ 
ally  returns  to  the  equilibrium  position,  then  its  motion  is  said  to  be 
damped.  If  insufficient  damping  is  available,  then  the  body's  motion 
will  become  divergent.  Divergent  oscillations  are  generally  undesire- 
able  and  usually  result  in  material  failure  of  a  mechanical  system. 
Damping  tends  to  take  energy  out  of  a  system.  Some  factors  which  cause 
damping  are  friction,  hydraulic  dampers,  springs,  etc.  Divergent 
motions  have  energy  added  to  the  dynamic  system.  An  example  of  this  is 
a  pilot-induced  oscillation  (PIO)  where  the  pilot's  control  movements 
are  in  the  same  direction  and  with  the  same  timing  as  the  aircraft's 
response. 

Neutral  stability  is  another  possible  behavior  of  a  disturoed  body. 
Here  the  body  remains  at  its  original  disturbed  state  or  oscillates  at  a 
constant  amplitude  about  the  equilibrium  state.  Static  stability  is 
necessary  for  dynamic  stability  to  exist,  but  the  converse  is  not  true. 
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0.  AXES  SYSTEMS 

To  establish  a  basis  for  the  discussion  of  unbalanced  forces  and 
moments  on  an  aircraft  and  that  aircraft's  subsequent  reaction,  a  ref¬ 
erence  coordinate  frame  must  be  established.  Figure  1-3  illustrates  the 
conventional  arrangement  of  perpendicular  axes  which  are  centered  at  the 
helicopter's  center  of  gravity.  The  directions  indicated  by  the  arrows 
are  the  positive  direction. 


Figure  1-3.  Axis  System  Notation 
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The  positive  directions  of  this  axis  system  for  the  X,  Y,  and  Z  axes 
are  forward,  right,  and  down,  respectively.  Thus  this  is  a  right-handed 
system.  (Note:  A  right-handed  rectangular  coordinate  system  derives 
its  name  from  the  analogy  that  a  right-threaded  screw  rotated  through  90 
degrees  from  OX  to  OY  will  advance  in  the  positive  Z  direction  (see 
Figure  1-4).  Forces  are  named  for  the  directions  along  which  they  act. 
Thus  an  X-force  is  one  acting  in  the  X-direction.  The  same  nomenclature 
system  applies  for  the  Y  and  Z  forces. 


Figure  1-4.  Right-Handed  Coordinate  System 


Rational  motion  also  occurs  about  the  X,  Y,  and  Z  axes.  These 
moments  are  termed  L,  M,  and  N.  I  is  the  rolling  moment  which  occurs 
about  the  longitudinal  (X)  axis.  A  roll  to  the  right  is  defined  as 
positive.  M  is  the  pitching  moment  which  occurs  about  the  lateral  (Y) 
axis.  A  nose-up  pitch  describes  a  positive  value  of  M.  The  yawing 
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N-moment  occurs  about  the  vertical  (Z)  axis.  A  positive  yaw  is  defined 
as  one  which  moves  the  nose  of  the  helicopter  to  the  right.  Table  1-1 
summarizes  the  axis  system  notation. 


TABLE  1-1.  Axis  System  Notation 


AXIS 

FORCE 

VELOCITY 

MOMENT 

ANGULAR 

MEASURE 

ANGULAR 

VELOCITY 

X 

X 

U 

_l 

0 

P 

Y 

Y 

V 

M 

0 

Q 

Z 

Z 

w 

N 

3 

R 

This  same  set  of  orthogonal  axes  can  be  referenced  in  various  ways, 
depending  on  the  needs  of  the  engineer.  Certain  particular  problems 
dealing  with  aircraft  stability  can  be  solved  more  easily  by  the  proper 
selection  of  axis  reference. 

Three  systems  of  axis  reference  are  generally  used: 

(1)  Gravity  Axis, 

(2)  Stability  Axis,  and 

(3)  Body  Axis. 

1.  Gravity  Axis 


In  this  system  the  Z-axis  of  the  helicopter  is  always  pointing 
at  the  center  of  the  earth,  and  the  X-axis  is  directed  along  the  horizon. 
The  gravity  axis  system  is  useful  for  linear  displacements  and  angular 


accelerations.  Certain  simplifications  in  the  stability  derivatives  can 
be  achieved;  however  when  helicopter  rotation  is  taken  into  account 
inertial  terms  and  products  of  inertia  must  have  lengthy  corrections. 

2.  Stability  Axis 

In  this  system  of  reference,  the  X-axis  is  aligned  with  the 
velocity  vector  and  is  positive  pointing  into  the  relative  wind.  The 
Z-axis  is  perpendicular  to  the  relative  wind  and  the  Y-axis  is  ortho¬ 
gonal  to  both,  forming  a  right-hand  system.  Using  the  stability  axis 
system  can  yield  great  simplifications  of  the  aerodynamic  terms.  This 
system  is  limited  to  small  disturbance  motions,  however,  because  the 
moment  of  inertia  terms  vary  and  thus  are  assumed  to  be  constant  in  the 
equations  of  motion. 

The  stability  axis  reference  system  has  the  useful  feature  of 
being  directly  applicable  to  wind  tunnel  results  which  are  commonly 
measured  parallel  with  and  perpendicular  to  the  wind. 

The  stability  axis  is  very  useful  for  fixed-wing  analysis, 
because  the  relative  wind  is  always  somewhat  directly  on  the  nose  of  the 
aircraft  and  varies  very  little  from  that  direction.  The  unique  capa¬ 
bility  of  a  helicopter  to  decelerate  from  forward  flight  to  hovering 
flight  (zero  forward  velocity)  means  that  the  X-  and  Z-axes  could  change 
by  as  much  as  90  degrees  as  the  relative  velocity  of  the  air  changes 
from  horizontal  (on  the  nose  of  the  helicopter)  to  vertical  (being 
pulled  down  through  the  rotor  disk).  Thus  the  usefulness  of  this  axis 
system  for  helicopters  is  limited  since  it  cannot  be  used  for  comparison 
purposes  over  the  whole  range  of  the  helicopter's  velocity. 


15 


3.  Body  Axis 

The  body  axis  system  aligns  the  X-axis  with  a  datum  line  on  the 
helicopter.  The  Z-axis  is  perpendicular  to  the  X-axis  and  is  directed 
out  the  bottom  of  the  aircraft  while  the  Y-axis  is  orthogonal  to  both. 
This  system  ensures  that  the  inertial  terms  in  the  equations  of  motion 
are  independent  of  the  flight  conditions. 

The  reference  system  of  body  axes  is  very  useful  when  studying 
helicopter  dynamics  because  velocities  and  accelerations  with  respect  to 
these  axes  are  the  same  as  those  that  would  be  measured  by  instruments 
in  the  helicopter  and  those  that  are  experienced  by  the  pilot. 
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II.  THE  EQUATIONS  OF  MOTION 

To  obtain  solutions  of  aircraft  stability,  some  quantitative  data 
must  be  made  available  to  the  engineer.  Naturally  some  formulae  or 
equations  would  be  helpful  when  trying  to  arrive  at  a  mathematical 
determination  of  the  stability  problem  at  hand. 

In  this  chapter  some  elementary  concepts  will  be  introduced  which 
will  lead  to  the  development  of  the  equations  of  motion  of  the  flight 
vehicle.  It  is  these  equations  of  motion  which  will  yield  the  numerical 
data  needed  for  problem  solutions. 

A.  LINEAR  MOTION 

Linear  motion  is  the  motion  of  an  object  along  a  line.  The  line  can 
either  be  straight  or  curved  so  that  linear  motion  can  be  further  sub¬ 
divided  in  rectilinear  (straight  line)  or  curvilinear  (along  a  curved 
line)  motion. 

1.  Rectilinear  Motion 

The  rectilinear  motion  of  a  particle  can  be  described  by  that 
particle's  position  on  a  straight  line  and  its  time  derivatives  of 
position.  To  quantify  this  motion  a  reference  point  must  be  selected. 
All  subsequent  measurements  of  the  particle's  motion  are  made  with 
respect  to  this  reference  point.  A  coordinate  system  well  suited  for 
recil inear  motion  is  the  rectangular  cartesian  coordinate  system. 

By  making  use  of  a  selected  coordinate  system,  p  marks  the 
generalized  coordinates  of  the  particle's  position.  The  distance,  s,  of 
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the  particle  from  the  selected  reference  point  is  the  difference  of  the 
coordinates  for  these  two  points.  If  the  origin  of  the  coordinate 
system  is  chosen  as  the  reference  point,  then  the  distance  to  the  parti¬ 
cle  simply  becomes  the  value  p.  If  some  point  other  than  the  origin  is 
chosen  as  the  reference  point,  then  the  distance  of  the  particle  from 
the  reference  point  is  p  -  r  where  r  is  the  coordinate  of  the  reference 
point  from  the  origin.  Should  the  particle  move  to  a  new  point  on  the 
line,  p 1  ,  it  would  then  be  at  a  different  distance  from  the  origin. 
This  difference  is  As  (see  Figure  2-1). 


Figure  2-1.  Measurement  Along  a  Line 


The  distance  is  given  as: 


As  =  p  -  p'  (2-1) 

The  rate  at  which  the  particle  travels  from  p  to  p‘  is  very  often  of 
interest.  Thus  the  average  velocity  of  the  particle  along  this  path  is 
the  distance  along  the  path  divided  by  the  time  it  took  to  cover  the 


di stance. 
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V  =  ££ 

avg  At 


(2-2) 


As  the  time  interval  becomes  very  small,  the  result  will  be  the  instant¬ 
aneous  velocity  of  the  particle  along  the  path. 


*  =  t'o  tt  (2'3) 

Both  average  and  instantaneous  velocities  will  be  used  in  subsequent 
calculations. 

Similarly,  it  can  be  shown  that  the  acceleration  of  the  particle 
is  the  rate  of  change  of  the  particle's  velocity. 


aavg 


A v 
At 


(2-4) 


In  the  same  manner  as  for  velocity,  the  instantaneous  acceler¬ 
ation  of  the  particle  can  be  shown  to  be: 


LIM  AV  _  dv  _  • 
3  "  t+0  At  dt  v 


(2-5) 


and  since  dv  =  ds/dt,  acceleration  can  also  be  expressed  as  the  second 
derivative  of  position  with  respect  to  time. 


_  ds/dt  _  d2s  _  .. 
a  dt  dt?  “  5 


(2-6) 


The  calculation  of  accelerations  become  important  because  they  lead 
directly  to  forces  on  the  airframe  via  Newton's  Second  Law: 

F  =  ma  (2-7) 
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2.  Curvilinear  Motion 

The  second  type  of  linear  motion  occurs  along  a  curved  path  and 
is  therefore  called  curvilinear  motion.  Plane  curvilinear  motion  occurs 
when  a  particle  that  is  moving  along  a  curved  path  remains  in  a  single 
plane.  For  this  case  the  position  in  which  the  particle  motion  occurs 
(see  Figure  2-2). 


Figure  2-2.  Curvilinear  Motion 


In  this  figure,  the  particle  has  moved  from  point  p  to  point  p1  along 
the  curved  path.  The  particle's  motion  takes  place  in  tne  plane  of  the 
paper  and  occurs  about  point  0.  The  position  of  the  new  point  p'  is 
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given  by  the  vector  addition  of  the  old  position  vector  at  p  (vector  r) 
plus  the  change  that  occurs  during  the  movement  (vector  Ar). 


As  was  the  case  with  rectilinear  motion,  the  average  velocity  of 
a  point  along  the  path  is  equal  to  the  time  rate  of  change  of  the 
point's  position  vector.  Thus 


v 

avg 


(2-8) 


If  the  time  interval  becomes  smaller  and  smaller,  the  Ar  vector  ap¬ 
proaches  tangency  to  the  curved  path.  In  the  limiting  case,  the 
velocity  will  approach  the  instantaneous  velocity  of  the  particle  along 
the  curved  path. 


_  LIM  Ar  _  dr 
t->-0  At  dt 


(2-9) 


If  the  speed  of  the  particle  increases  as  it  moves  from  p  to  p' , 
then  the  particle  is  accelerating.  The  same  steps  that  were  used  to 
develop  the  instantaneous  velocity  can  be  followed  to  find  the  acceler¬ 
ation  term.  The  result  is 


_  LIM  Av  _  dv  _  df_r 
t-*-0  At  ~  at  dt? 


B.  ANGULAR  MOTION 


(2-10) 


While  linear  motion  can  be  used  to  describe  a  significant  number  of 
motions  that  are  commonly  encountered  in  physics,  motion  can  also  be 
described  in  another  way.  Angular  measurements  to  a  moving  point  from  a 
fixed  point  are  commonly  used  when  dealing  with  rotating  systems. 
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The  description  of  angular  motion  will  be  limited  to  that  occuring 
in  a  plane.  Before  this  angular  motion  can  be  described,  some  reference 
point  must  be  selected.  A  point  in  the  plane  of  motion  is  a  good  refer¬ 
ence  point  and  an  axis  perpendicular  to  the  plane  of  motion  serves  as 
the  axis  about  which  the  angular  motion  occurs.  This  becomes  the  axis 
of  rotation. 

Before  angular  measurement  can  begin,  a  reference  axis  in  the  plane 
of  rotation  must  be  chosen.  All  angular  positions  will  be  measured  with 
respect  to  this  reference  axis.  The  reference  axis  itself  is  chosen 
arbitrarily  and  it  does  not  rotate.  Although  the  system  being  measured 
may  rotate,  the  reference  axis  remains  fixed. 

Angular  displacements  are  measured  in  degrees  or  radians  using  the 
symbol  8.  By  convention,  a  positive  angular  displacement  is  in  a 
counterclockwise  direction  from  the  reference  axis  (see  Figure  2-3). 


0 


Figure  2-3.  Angular  Measure 
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According  to  the  foregoing  rules,  the  angular  measurement  of  line  OB 
is  at  some  angle  9  from  the  arbitrarily  chosen  reference  axis,  OA.  The 
axis  of  rotation  is  perpendicular  to  the  paper  at  point  0.  The  angular 
rate  at  which  the  line  moves  from  position  OA  to  OB  in  Figure  2-3  is  the 
time  rate  of  change  of  its  angular  position.  As  was  the  case  with 
previous  descriptions  of  linear  velocity,  the  angular  velocity,  w,  is 
the  first  time  derivative  of  the  angular  position. 

“  =  f  *  ®  (Z-U) 

Similarly,  angular  acceleration,  a ,  is  the  time  derivative  on  angular 
velocity  or  the  second  time  derivative  of  angular  position. 

<*  =  at  =  *  =  0  (2-12) 

A  common  case  of  motion  combines  both  curvilinear  and  angular 
motions.  This  case  of  motion  occurs  when  a  particle  moves  around  a 
fixed  point  at  a  fixed  distance  from  that  point.  An  example  of  this  is 
the  arc  traced  out  by  the  helicopter's  rotor  blade  tips. 

As  a  helicopter  blade  sweeps  the  air,  all  points  on  the  blade  are 
rotating  at  the  same  angular  velocity.  However  each  point  at  a  differ¬ 
ent  radius  from  the  center  of  the  axis  of  rotation  has  a  different 
curvilinear  velocity  or  speed.  This  can  be  readily  understood  by  ex¬ 
amining  Figure  2-4.  and  R2  are  two  points  located  on  the  same 
rotating  line  (for  example,  a  rotor  blade).  During  a  single  rotation  of 
the  rotor  disk,  R2  travels  a  greater  distance  than  Rj.  The  distance 
traveled  by  each  point  is  the  circumference  of  the  circle  traced  out  by 
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Figure  2-4.  Dependence  of  Linear  Velocity  on  Radius 


the  point.  The  path  traced  out  by  R2  is  longer  than  that  traced  out  by 
Rt.  Since  both  points  complete  one  revolution  of  the  circle  in  the  same 
amount  of  time,  R2  must  travel  faster  than  Rx. 

The  velocity  of  a  particle  rotating  about  a  fixed  point  is  given  by 
both  its  angular  velocity  and  distance  from  the  axis  of  rotation  and  is 
developed  in  the  following  manner  (refer  to  Figure  2-5).  The  instant¬ 
aneous  velocity  of  the  rotating  point  is  the  limit  as  time  approaches 
zero  of  the  change  in  the  particle's  position  vector  divided  by  tne 
average  time  over  which  this  position  change  occurs. 

_  1  im  Sr 
t-0  4t 
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v 


(2-13) 


Figure  2-5.  Linear  Velocity  of  a  Rotating  Point 


where  Ar  =  r  times  the  change  in  angular  displacement.  Thus 


v 


1  i  m  rA6 
t->0  At 


(2-14) 


r  is  measured  in  feet  or  meters  and  w  is  measured  in  radians  per  second. 

It  can  be  seen  from  equation  (2-14)  that  a  particle  located  at  half 
the  distance  from  the  rotor  hub  to  the  tip  of  the  rotor  is  moving  at 
half  the  linear  velocity  of  a  particle  located  at  the  rotor  tip. 


C.  FORCES  AND  MOMENTS 


1.  Forces 

No  discussion  of  forces  would  be  complete  without  giving  credit 

to  Sir  Issac  Newton  and  his  statements  of  the  basic  laws  governing  the 

motion  of  a  particle  [Ref.  2].  They  are: 

Law  I.  A  particle  remains  at  rest  or  continues  to  move  in 
a  straight  line  with  a  uniform  velocity  if  there  is  no 
unbalanced  force  acting  on  it. 

Law  II.  The  acceleration  of  a  particle  is  proportional  to 
the  resultant  force  acting  on  it  and  is  in  the  direction  of 
the  force. 

Law  III.  The  forces  of  action  and  reaction  between  inter¬ 
acting  bodies  are  equal  in  magnitude,  opposite  in  direction, 
and  col  linear. 

The  word  force  is  mentioned  in  all  three  of  Newton's  laws. 
Therefore  it  would  seem  that  force  plays  a  large  part  in  the  movement  of 
a  particle  (or  of  an  object,  if  tne  center  of  gravity  of  that  object  is 
considered  to  be  a  particle).  Force  can  be  defined  as  the  action  of  one 
body  on  another.  A  force  acting  in  a  direction  tends  to  move  an  object 
it  acts  upon  on  the  same  direction,  according  to  Law  II. 

It  should  be  noted  that  forces  are  vector  quantities.  That  is, 
they  are  composed  of  a  magnitude  and  a  direction.  The  direction  of  a 
force  becomes  very  important  when  attempting  to  predict  the  reactions  of 
objects  to  forces. 

The  importance  of  the  vectorial  nature  of  forces  cannot  be  over¬ 
emphasized.  For  example,  what  is  the  resultant  of  a  50-pound  force 
applied  to  another  50-pound  force?  The  correct  answer  could  be  zero 
pounds,  100  pounds,  or  anywhere  in  between  depending  on  whether  the 
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forces  were  directly  opposite  to  each  other,  in  the  same  direction,  or 
at  some  other  angular  position  between  these  two  extremes. 

Newton's  second  law  is  of  primary  interest  for  the  moment. 
Another  way  of  stating  the  second  law  is  that  the  resultant  force  acting 
on  a  particle  is  proportional  to  the  time  rate  of  change  of  the  momentum 
of  the  particle  and  that  this  change  is  in  the  direction  of  the  force. 
The  proportionality  factor  is  the  mass  of  the  particle  in  question. 

Both  statements  of  Newton's  second  law  lead  to  the  same  result 
in  equation  form: 


L?  _  dmv  _  mdv  _  -  /o-ic-v 
F  -  “dt  -  -at  •  ma  (2  16) 

It  will  be  assumed  that  the  mass  of  the  particle  does  not  change 
during  the  time  interval  dt.  Thus  dmv/dt  can  be  written  mdv/dt.  Also 
note  that  I  symbol  was  used.  The  resultant  acceleration  of  the  particle 
is  equal  to  the  resultant  vector  sum  of  the  ^orces  acting  on  the  parti¬ 
cle.  Thus  many  different  forces  acting  on  an  object  can  be  vectorial ly 
resolved  into  one  force  (see  Figure  2-6). 


Figure  2-6.  Vector  Summation 
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The  particle  will  -espond  in  the  same  way  no  matter  if  we  deal 
with  the  several  individual  forces  or  with  their  resultant,  but  the 
picture  is  greatly  simplified  when  dealing  with  only  a  single  resultant 


fc 


force. 

Both  the  vectorial  nature  of  forces  and  the  capability  to  sum 
individual  forces  can  be  used  with  the  axis  systems  previously  described 
for  the  helicopter.  Being  a  complex  machine,  the  helicopter  can  have 
many  different  forces  from  different  sources  acting  at  various  points  on 
the  body  of  the  aircraft.  These  forces  can  be  broken  down  into  compon¬ 
ents  along  the  aircraft  axes  and  summed  such  that  there  now  exists  three 
mutually  perpendicular  forces  acting  at  the  center  of  gravity  of  the 
aircraft.  The  problem  of  the  helicopter's  reaction  to  these  various 
forces  is  now  greatly  simplified  (see  Figure  2-7). 

X 

A 
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Figure  2-7.  Components  of  a  Vector  in  a  Plane 


Equation  (2-15)  can  be  written  for  each  axis  direction: 


Efx  =  max  ;  =  maY  ;  =  maZ  (2-16) 

As  an  example,  a  single  force  not  lying  in  any  single  plane  would  have 
components  in  the  X,  Y,  and  Z  directions.  The  components  of  this  force 
would  be  determined  and  would  be  respectively  included  in  each  of  the 
three  force  equations  above.  If  a  force  acts  entirely  along  a  single 
axis,  the  Z-axis  for  example,  it  would  have  no  components  in  the  other 
two  directions,  X  or  Y  in  this  case. 

Units  of  measure  for  force  are  listed  below  in  Table  2-1. 


Table  2-1.  Units  of  Measure  for  Force 


SYMBOL 

ENGLISH 

ENGINEERING 

SYSTEM 

SI 

SYSTEM 

FORCE 

F 

Pounds  Force 

Newtons 

MASS 

m 

Slugs 

Kilograms 

ACCELERATION 

a 

ft/ sec2 

m/sec2 

Since  the  consideration  of  variable  masses  will  not  be  con¬ 
sidered  (the  mass  of  a  helicopter  will  not  change  over  the  time  periods 
under  consideration  here)  it  can  be  seen  that  the  acceleration  of  a  body 
is  directly  proportional  to  the  force  acting  on  the  body. 
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2.  Moments 


Forces  applied  to  a  particle  create  only  linear  motion  of  the 
particle.  This  is  so  because  a  particle  is  considered  to  be  a  very 
small  quantity  and  any  forces  applied  must  necessarily  act  at  its  center 
of  mass.  The  phenomenon  of  angular  acceleration  occurring  simultaneous¬ 
ly  with  linear  acceleration  arises  when  a  force  is  applied  to  an  object 
(something  larger  than  a  particle).  If  the  direction  vector  of  a  force 
applied  to  an  object  does  not  pass  directly  through  the  object's  center 
of  mass  then  the  object  will  start  to  rotate  about  its  center  of  mass 
because  a  moment  has  been  created. 

A  moment  is  a  force  applied  at  a  distance  from  an  axis  around 
which  an  object  rotates.  If  the  object  is  located  in  free-space,  it 
will  rotate  about  its  center  of  mass.  If  the  object  is  pinned  somehow, 
the  pin  will  act  as  the  pivot  point  about  which  the  object  will  rotate. 

In  equation  form,  for  a  fixed- axis  system: 

M  =  Fr  (2-17) 

where  M  is  the  moment  in  foot-pounds,  F  is  the  force,  and  r  is  the  per¬ 
pendicular  distance  from  the  object's  center  of  mass  to  the  force  vector 
(see  Figure  2-8). 

Just  like  a  force,  a  moment  is  a  vectorial  quantity  but  it  is  a 
rotational  vector  in  this  case.  As  with  angular  velocity,  a  positive 
moment  is  defined  in  the  counterclockwise  direction. 

In  the  evaluation  of  the  effects  external  moments  have  on  the 
•notion  of  the  helicopter,  the  moment  of  momentum,  H.  is  considered. 
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Figure  2-8.  Generation  of  a  Moment 


By  definition,  the  moment  of  momentum  of  a  small  part  on  the  aircraft, 
dm,  is: 

6h  =  r  x  V  dm  (2-18) 

Taking  the  derivative  of  the  above  yields 

^  5h  =  (r  x  V)  dm  =  ^  x  V  dm  +  r  x  ^  6m  (2-19) 

Considering  a  small  chunk  of  mass  on  a  rotating  body,  the  velo¬ 
city  of  this  piece  of  mass  is: 

V  =  V  +  ^  (2-20) 

c  dt 

Where  dr/dt  is  the  rotational  velocity  of  the  body  and  is  the  velo¬ 
city  of  the  object's  center  of  mass. 


A1  so, 


~r  x  ai 5m  = 


r  x  6F  =  6G 


(2-21) 


This  is  true  since  dV/dt  6m  =  6F  is  just  F  =  ma,  while  6G  is  the  moment 
about  the  center  of  mass  of  the  object  produced  by  force  5F. 

Substituting  equations  (2-20)  and  (2-21)  into  equation  (2-19) 

gives: 

5G  =  (6h)  -  (V  -  V  )  x  V  6m  (2-22) 

where  G  is  the  linear  momentum  of  a  particle,  mv. 

Since  V  x  V  =  0,  this  becomes 


6G  =  ^  (5h)  +  Vc  x  V  dm 


(2-23) 


If  equation  (2-23)  is  summed  for  all  small  mass  elements  on  the 
helicopter  the  resultant  mass  is  the  total  mass  of  the  helicopter  and 
the  resultant  velocity  is  V  .  As  shown  before,  V  x  V  =  0  and  equation 
(2-23)  reduces  to 


s  _  dh 
G  “  dt 


(2-24) 


This  equation  states  that  the  angular  momentum  of  an  object  is  changed 
when  a  moment  is  applied  to  the  object. 
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D.  CORIOLIS  FORCES  AND  MOMENTS 


Coriolis  forces  arise  from  the  acceleration  produced  when  a  particle 
moves  along  a  path  in  a  plane  which  is  itself  rotating.  A  flapping 
blade  is  subject  to  a  coriolis  acceleration  in  the  plane  of  rotation. 

When  a  rctor  blade  flaps,  its  moment  of  inertia  about  the  rotational 
axis  changes.  This  can  be  seen  by  considering  the  center  of  mass  of  the 
blade  being  rotated  to  an  increased  flapping  angle,  j$  +  AB.  As  the 
blade  flaps  up,  the  center  of  mass  moves  closer  to  the  rotational  axis 
(see  Figure  2-9'. 


HoriZcW 


Figure  2-9.  Motion  of  Flapping  Blade 


From  the  law  of  the  conservation  of  angular  momentum,  the  blade  will 
experience  an  accelerating  force  if  the  center  of  mass  of  the  blade 
moves  closer  to  the  rotational  axis  (flaps  up  from  A  to  B)  or  a  retard¬ 
ing  force  if  the  center  of  mass  moves  farther  away  (flaps  down  from  B  to 
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A).  These  forces  are  manifested  as  vibrations  at  the  blade  root,  and 
for  every  harmonic  of  flapping  there  is  an  appropriate  inplane  Coriolis 
vibration. 

To  find  the  coriolis  force  the  following  steps  are  used  [Ref.  3]. 
Point  B  has  two  components  of  velocity  which  are  of  interest,  dr/dt 
along  the  blade  projection  (inwards)  and  uir  cos  (p  +  Ap)  (in  the  plane 
of  rotation).  Taking  components  of  velocity  on  the  line  perpendicular 
to  the  blade  at  P2,  (see  Figure  2-10)  results  in: 

AS  +  Qr  cos  (p  +  Ap)  (2-24) 

and  the  difference  in  velocity  perpendicular  to  the  original  projection 
between  P2  and  Px  is: 

^  A8  +  Qr  (cos  p  +  AP)  -  Qr  cos  p 

(2-25) 

=  {  ^  (r  cos  p)f A0  +  Qr  (cos  p  cos  AP  -  sin  p  sin  Ap)  -  Qr  cos  p 


Figure  2-10.  Flapping  Blade  in  Plane  of  Rotation 
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Equation  (2-25)  was  arrived  at  by  using  the  trigonometric  identity 

cos  0  +  Ap  =  cos  p  cos  Ap  -  sin  p  sin  Ap  (2-26) 

Using  small  angle  theory  cos  p  1  and  sin  p  p.  This  will  result 
in 

AV  =  r  sin  p  jj|  A8  -  Or  p  Ap  (2-27) 


Dividing  through  by  At  to  find  the  acceleration  yields: 


AV  _ 
At  " 


r«  A9 
rp  dt  At 


Qr  p 


At 


and  since 


$  =  Q  and  $  = 
At  At 


(2-28) 


acor  =  ‘  <"£  PG  ‘  Qr  pp  (2-29) 

which  gives  the  final  form  of  the  Coriolis  acceleration: 

acor  =  ’  2  PP  (2-30) 

If  the  Coriolis  force  is  desired,  simply  multiply  the  mass  of  the 
blade  by  the  acceleration  obtained  above. 

Fcor  *  '  2  "b  ^  «  <2-31) 

This  Coriolis  force  can  be  considered  to  act  at  the  center  of  the 
blade  and  it  produces  a  moment  about  the  hub  of  the  rotor.  The  moment 
produced  therefore  is  the  force  times  the  distance  through  which  the 
force  acts, 
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r 


(2-32) 


f-1  „  =  2  M.  r2  Q  pp 
cor  b 

A  useful  relation  can  be  used  to  shorten  this  equation.  The  moment 
of  inertia  of  an  object  is  defined  as 

I  =  y'y2dm  (2-33) 

where  y  is  the  distance  from  an  axis  and  dm  is  an  elementary  particle  of 
mass.  If  the  total  mass  of  the  object  and  its  center  of  mass  location 
is  known,  the  moment  of  inertia  can  be  expressed  as 


I  =  k2m  (2-34) 

where  k  is  the  distance  from  the  axis  to  the  center  of  mass  and  m  is  the 

total  mass.  Thus,  substituting  r  for  k  in  equation  (2-34)  and  substi¬ 

tuting  equation  (2-34)  into  equation  (2- 32)  yields 

Mcor  =  2  I  0  PP  (2-35) 

NOTE:  The  Coriolis  force  just  derived  is  not  the  only  Coriolis  force 

that  arises  from  the  dynamic  notions  of  the  helicopter.  It  is 
perhaps  the  most  significant  and  easily  understood  Coriolis 
force,  but  it  must  be  pointed  out  that  many  Coriolis  forces  must 
be  accounted  for  when  a  very  rigorous  analysis  of  the  heli¬ 
copter's  motions  is  conducted.  Other  Coriolis  forces  will  arise 
from  aeroelastic  effects  of  the  rotor  blades,  differences  in 
motion  between  the  helicopter  rotor  and  fuselage  and  also  from 
fuselage  aeroelastic  effects.  Many  of  these  latter  Coriolis 
forces  will  be  very  small  and  will  not  be  apparent  when  con¬ 
sidering  the  helicopter  as  a  rigid  body. 
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E.  RIGID  BODY  EQUATIONS 

Now  that  some  of  the  background  concerning  the  nature  of  forces  and 
moments  has  been  presented,  it  remains  to  put  these  ideas  to  use  for  the 
solution  of  problems  dealing  with  helicopter  stability.  Simply  put,  the 
various  forces  and  moments  that  act  on  the  helicopter  are  resolved  into 
components  which  act  along  or  about  one  of  the  principal  axes  of  the 
helicopter.  The  various  forces  and  moments  to  be  considered  are  the 
factors  which  give  rise  to  the  vehicle's  motion. 

According  to  Reference  4,  certain  assumptions  concerning  the  heli¬ 
copter's  motions  and  references  will  have  to  be  made  in  order  to 
establish  the  ground  rules  for  further  analysis  and  also  to  simplify  the 
systems  of  equations  that  will  arise. 

Assumption  1:  The  helicopter  is  a  rigid  body. 

A  rigid  body  is  one  in  which  motions  between  individual  mass  ele¬ 
ments  that  make  up  the  body  do  not  occur.  In  this  way,  distances 
between  specified  joints  in  the  helicopter's  fuselage  remain  fixed. 
Thus  no  bending  or  twisting  of  the  fuselage  is  considered  as  the  heli¬ 
copter  moves  through  space. 

This  assumption  allows  the  helicopter's  motion  through  space  to  be 
described  as  the  rectilinear  motion  of  the  center  of  mass  of  the  air¬ 
craft  and  by  the  curvilinear  motion  of  the  same  point.  In  reality,  the 
fuselage  does  bend  and  twist  to  a  certain  degree  in  flight.  The  dy¬ 
namics  of  the  structure  are  changed  when  these  aeroelastic  effects  are 
considered.  The  aeroelastic  effects  serve  to  greatly  increase  the 
degrees  of  freedom  that  must  be  considered  when  analyzing  the  equations 
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of  motion.  A  rigid  body  will  be  assumed  here  since  aeroelasticity  is 
beyond  the  purpose  of  this  paper  and  consequently  the  solution  to  gener¬ 
ated  problems  will  not  be  overly  complicated. 

Assumption  2:  The  earth  is  considered  to  be  flat  and  is  fixed  in 
space. 

This  assumption  is  made  to  negate  any  minor  corrections  that  may 
otherwise  have  to  be  made  for  the  gravity  constant  or  for  a  moving  in¬ 
ertial  reference  frame.  Negligible  error  is  introduced  by  this  assump¬ 
tion  since  the  altitudes,  air  speeds,  and  time  lengths  under  consider¬ 
ation  are  small. 

Assumption  3:  The  helicopter  is  assumed  to  be  in  a  trimmed  level 
flight  condition,  or  in  a  hover.  Small  perturbations 
in  the  helicopter's  motion  are  then  considered. 

This  assumption  allows  the  linearizing  of  normally  nonlinear  re¬ 
sponses  by  considering  only  small  increments  of  motion.  Simplification 
of  the  required  equations  used  to  solve  for  the  helicopter’s  motion  is 
the  result.  While  the  analyses  will  be  performed  for  very  small  changes 
in  motion,  the  results  can  be  extended  to  larger  motions  without  very 
much  loss  in  accuracy. 

Other  mathematical  benefits  arising  from  this  assumption  include  the 
small  angle  approximations  for  sine  and  cosine: 

(1)  sina  =  a 

(2)  cosa  =  1 

where  a  is  measured  in  radians,  and  the  fact  that  the  power  of  any  small 
increment  can  be  ignored  when  compared  to  the  original:  e2  =:  0  when  e  is 
small. 
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Assumption  4:  The  helicopter  is  trimmed  in  steady,  level  flight  and 
the  longitudinal  forces  and  moments  due  to  lateral 
perturbations  from  trim  are  negligible. 

This  assumption  discounts  any  initial  angular  velocities  in  roll  or 
yaw  (P0,  R0)  and  any  initial  lateral  velocity  (V0).  Furthermore,  any 
lateral  perturbations  that  do  arise  will  not  affect  the  longitudinal 
response  of  the  aircraft.  This  last  statement  is  important  since  the 
longitudinal  and  lateral  equations  of  motion  are  thereby  decoupled. 

Assumption  5:  The  X-Z  plane  is  a  plane  of  symmetry. 

This  is  usually  the  case  for  most  flight  vehicles.  Although  this  is 
not  true  for  many  helicopters,  very  little  error  is  introduced  into  the 
final  equations  of  motion  when  this  assumption  is  made.  With  symmetry 
and  with  Assumption  4  considered  all  rolling  moments,  yawing  moments, 
and  side  forces  are  reduced  to  zero.  Consequently  the  longitudinal 
equations  of  motion  can  be  described  by  just  three  equations:  those 
dealing  with  longitudinal  and  vertical  forces  and  with  the  pitching 
moment. 

Combining  these  five  assumptions  makes  it  possible  to  describe  the 
helicopter's  motion  with  just  six  equations.  This  is  so  because  the 
assumptions  allowed  many  additional  degrees  of  freedom  to  be  dropped 
from  consideration.  Consequently,  three  force  equations  and  three 
moment  equations  are  used  to  describe  the  helicopter's  motion. 

Note  in  particular  Assumption  4.  Since  the  longitudinal  ana  lateral- 
directional  motions  are  now  decoupled,  only  those  forces  and  moments 
which  affect  the  longitudinal  motions  of  the  helicopter  need  be  con¬ 
sidered.  These  turn  out  to  be  the  X  and  Z  forces  and  the  M  moment. 
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1.  Three  Force  Equations 

Since  forces  act  in  straight  lines,  it  is  very  convenient  to 
resolve  the  components  of  given  forces  acting  on  the  helicopter  along 
the  X,  Y,  and  Z  axes.  Only  external  forces  on  the  aircraft  will  be 
considered.  Internal  forces  must  necessarily  be  opposed  by  other  in¬ 
ternal  forces  such  that  the  sum  total  of  all  the  internal  forces  on  the 
helicopter's  fuselage  is  zero.  While  internal  forces  may  play  an  im¬ 
portant  part  in  airframe  structural  considerations,  they  play  no  part  in 
the  analysis  of  the  helicopter's  motions.  The  external  forces  that  will 
be  considered  arise  from  three  principal  sources:  gravity  forces,  in¬ 
ertia  forces,  and  aerodynamic  forces. 

2 .  Three  Moment  Equations 

The  same  principles  that  were  used  for  the  determination  of  the 
three  force  equations  are  also  used  for  moments.  Here,  however,  instead 
of  the  X,  Y,  and  Z  axes  acting  as  directions  along  which  forces  a.  e 
measured,  the  axes  act  as  lines  around  which  the  moments  turn.  No¬ 
menclature  of  the  moments  will  be  briefly  reviewed: 

L-Moment:  Rotates  about  the  X-axis;  positive  is  right  wing  down. 

M- Moment:  Rotates  about  the  Y-axis;  positive  is  nose  up. 

N-Moment:  Rotates  about  the  Z-axis;  positive  is  nose  right. 

F.  EXPANSION  OF  FORES  AND  MOMENTS 

The  equatons  for  the  generalized  aerodynamic  forces  and  moments 
acting  on  an  aircraft  have  commonly  been  written  in  the  form 

A  =  CA  'i  p  V2S  (2-36) 
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where  A  stands  for  any  desired  force  or  moment,  CA  is  a  nondimensional 
coefficient,  p  is  air  density,  V  is  the  steady-state  velocity  of  the 
vehicle,  and  S  is  the  surface  area.  It  can  be  seen  by  this  equation 
that  the  aerodynamic  forces  and  moments  generated  on  the  helicopter  are 
dependent  upon  the  density  of  the  air  through  which  the  vehicle  is 
flying  and  also  the  velocity  of  the  aircraft  relative  to  the  air. 

The  nondimensional  coefficients,  C^,  are  also  dependent  on  Reynolds 
and  Mach  numbers,  angles  of  attack  and  sideslip,  and  linear  and  angular 
velocity  and  their  derivatives.  If  the  aerodynamic  forces  are  consider¬ 
ed  to  be  continuous  functions  of  all  these  variables,  each  of  the  forces 
and  moments  (X,  Y,  Z  and  L,  M,  N  respectively)  can  be  expressed  in  terms 
of  the  variables  by  expanding  the  terms  in  a  Taylor  series  [Ref.  4], 

A  Taylor  series  for  the  effect  of  forward  velocity  changes  has  the 
form: 


F 
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(2-37) 


All  partial  derivative  terms  starting  with  the  second  derivative  ar'e 
higher-order  terms  and  can  be  neglected  by  means  of  Assumption  3  in  the 
previous  section  without  changing  the  accuracy  of  the  solution.  This 
will  help  to  simplify  the  equation  by  eliminating  many  terms  which  serve 
to  complicate  the  problem  but  contribute  very  little  to  the  value  of  the 
final  solution. 

The  force  expression  is  now  reduced  to  one  showing  the  initial  trim 
condition  and  all  changes  to  the  force  resulting  from  the  first  partial 
derivative  of  force  with  respect  to  major  variables  and  their  deriva¬ 
tives  when  applicable. 
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(2-38) 


F  =  Fo  ♦  §3  4U  *  I  “  *  I?  ie 
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Since  the  change  in  force  from  a  steady  flight  condition  due  to 
perturbations  is  desired  and  not  the  total  force,  the  trim  force  can  be 
subtracted  from  both  sides  of  the  equation.  The  result  is: 


4F  *  I  4u  *  ti  42  *  §5  49  *  ^  “  *  fi  ai  * 


(2-39) 
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III.  CALCULATION  OF  THE  STABILITY  DERIVATIVES 


Probably  the  most  limiting  factor  of  expanding  the  helicopter's 
utility  to  a  multi-mission  role  has  been  its  relatively  slow  cruise 
speed.  Therefore  much  emphasis  has  been  directed  in  recent  years  to¬ 
wards  increasing  the  maximum  forward  flight  velocity  of  helicopters. 
The  stability  response  of  helicopters  in  hover  and  in  slow  flight  can 
generally  be  classified  as  one  consisting  of  two  heavily  damped  roots 
and  a  divergent  oscillatory  motion.  These  characteristics  can  (and 
usually  do)  change  as  the  aircraft  increases  its  forward  speed. 

Even  though  the  response  of  a  helicopter  in  a  hover  is  generally 
known,  its  response  at  high  forward  flight  velocities  cannot  accurately 
be  predicted.  Therefore,  the  engineer  must  be  able  to  calculate  the 
stability  derivatives  for  the  aircraft  in  question  and  at  the  proper 
flight  condition  (altitude  and  airspeed)  so  that  he  may  solve  for  the 
helicopter's  modes  of  response.  The  equations  and  procedures  presented 
in  this  chapter  outline  how  the  calculation  of  those  stability  deriva¬ 
tives  are  carried  out.  The  charts  used  for  this  procedure  apply  to  the 
regime  of  high  forward  flight  velocity  with  an  advance  ratio  of  0.3  or 
greater. 

A.  TYPES  AND  USES 

Both  dimensional  and  non-dimensional  stability  derivatives  are  used 
for  the  solution  of  problems  dealing  with  aircraft  stability  and  control. 


Each  type  derivative  has  certain  useful  properties.  By  using  non-dimen¬ 
sional  derivatives,  the  stability  characteristics  of  aircraft  can  be 
compared  regardless  of  size. 

Non-dimensional  derivatives  are  concerned  with  force  and  moment  co¬ 
efficients  and  with  non-dimensional  velocity  and  time.  The  real  advan¬ 
tage  to  using  non-dimensional  stability  derivatives  occurs  when 
comparing  the  stability  values  between  different  sizes  of  the  same 
aircraft,  for  example  comparing  data  obtained  from  a  one-tenth  scale 
model  in  a  wind  tunnel  with  data  expected  or  obtained  from  flight  tests 
of  the  full  scale  aircraft.  This  capability  of  non-dimensional  deriva¬ 
tives  to  allow  the  comparison  of  data  in  this  manner  makes  possible  the 
prediction  of  stability  characteristics  of  aircraft  based  on  tests 
performed  on  scale  models.  It  is  immediately  apparent  that  this  is  an 
economical  benefit  when  the  al ternati ve  ’  i s  to  perform  all  tests  and 
development  on  full  scale  aircraft. 

Dimensional  derivatives  are  good  for  measuring  direct  forces  and 
moments  of  the  aircraft.  The  use  of  dimensional  stability  derivatives 
lead  directly  to  numerical  coefficients  in  the  sets  of  simultaneous 
differential  equations  describing  the  real  time  dynamics  of  the  airframe. 
By  analyzing  the  dynamic  response  numbers  thus  obtained,  the  heli¬ 
copter's  stability  characteristics  can  be  ascertained.  Once  this  is 
known,  the  amount  of  control  necessary  to  reach  the  desired  flying 
qualities  can  be  added. 

Care  should  be  taken  when  comparing  the  values  of  stability  deriva¬ 
tives  received  from  different  sources.  It  is  important  to  note  that 
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both  dimensional  and  non-dimensional  derivatives  are  widely  used 
throughout  the  industry  and  that  different  methods  may  be  used  to  non- 
dimensional  ize  or  normalize  the  derivatives  for  a  particular  aircraft. 

Thus  if  the  values  of  the  stability  derivatives  for  a  given  aircraft 
from  one  source  are  found  to  differ  from  those  values  received  from 
another  source,  the  most  likely  cause  for  the  difference  is  that  a 
different  system  of  non-dimensionalizing  or  normalizing  the  derivatives 
has  been  used  by  each  source. 

It  should  be  noted  that  one  common  method  used  for  dimensional 
derivatives  is  to  normalize  the  force  derivatives  by  the  mass  of  the 
aircraft  and  the  moment  derivatives  by  the  aircraft's  moment  of  inertia 
about  the  pitching  axis.  This  method  is  used  here. 

In  this  chapter  the  procedures  will  be  outlined  for  calculating  the 
stability  derivatives  Xu,  X^,  Zu,  Zw,  Mu,  and  M  .  These  derivatives 
are  needed  to  solve  the  stability  determinant  which  will  yield  the  modes 
of  motion  of  the  helicopter.  These  procedures  are  presented  in  more 
detail  in  [Ref.  5]. 

Initially,  certain  geometric  data  must  be  known  about  the  helicopter. 
The  forward  flight  velocity  at  which  the  stability  condition  will  be 
evaluated  must  also  be  chosen,  and  the  trim  condition  of  the  helicopter 
computed.  Since  the  aerodynamic  forces  which  act  on  the  helicopter  are 
dependent  on  flight  velocity,  the  trim  condition  and  subsequent  stabil¬ 
ity  derivatives  will  be  different  for  each  velocity  evaluated. 

After  the  trim  conditions  of  the  helicopter  are  found  at  the  desired 
airspeed,  the  isolated  derivatives  are  determined  and  corrected  for  main 
rotor  solidity  to  yield  the  local  derivatives.  From  these,  the  total 
derivatives  may  be  calculated. 


45 


8.  DEFINITION  OF  TERMS 


The  following  definitions  will  be  used: 

a  =  Lift  curve  slope  of  the  rotor  blade 

a0  =  Blade  coning  angle  (radians) 

at  =  Longitudinal  flapping  angle  (radians) 

AR  =  Aspect  ratio 

Av  =  Fuselage  frontal  area  (sq  ft) 

*FUS 

A,  =  Fuselage  planform  area  (sq  ft) 

^FUS 

b  =  Number  of  rotor  blades 

bt  =  Lateral  flapping  angle  (radians) 

c  -  Blade  chord  (ft) 

CD  =  Profile  drag  coefficient 

CQ  =  Drag  coefficient  =  DApV2$ 

Cg'  =  Drag  coefficient  of  the  main  rotor  =  D/T.F. 

=  Lit  '  coefficient  =  LApV2S 

C^'  =  Lift  coefficient  of  the  main  rotor  =  L/T.F. 

Cg  =  Rotor  torque  coefficient  =  Q/T.F.  R 

D  =  Aerodynamic  drag  force  (lbs) 

e  =  Blade  hinge  offset 

Ib  =  Blade  moment  of  inertia  about  flapping  hinge  (slug-sq  ft) 

Iy  =  Aircraft  moment  of  inertia  about  the  Ditching  axis 

K  =  Dcwnwash  interference  factor 
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L  =  Aerodynamic  lift  force  (lbs) 

1^  =  Longitudinal  moment  arm,  positive  when  the  point  of  appli¬ 

cation  of  the  force  ''ector  is  forward  from  the  C.G.  position 
(ft) 

U  =  Normal  (vertical)  moment  arm,  positive  when  the  point  of 
application  of  the  force  vector  is  below  the  C.G.  position 
(ft) 

M  =  Pitching  moment  of  an  aircraft  component  (ft- lb) 

Mu,  Mw>  Mq  =  Pitching  moment  total  derivatives 

Mg  =  First  moment  of  blade  mass  about  the  flapping  hinge  (slug-ft) 
Mj  =  Mach  number  of  advancing  blade  tip 

Q  =  Rotor  torque  (ft- lbs) 

q0  =  Dynamic  pressure  =  ^pV2  (Ib/sq  ft) 

R  =  Rotor  radius  (ft) 

S  =  Area  of  an  aerodynamic  surface  (sq  ft) 

T  =  Rotor  thrust  (lbs) 

T. F.  =  Thrust  factor  =  pnR2(QRj2  (lbs) 

V0  =  Steady  state  or  trim  value  of  velocity  (ft/sec) 

V<-  =  Velocity  of  sound  in  standard  atmosphere  (ft/sec) 

W  =  Aircraft  gross  weight  (lbs) 

X  =  Longitudinal  force  along  the  body  X-axis  (lbs) 

X^,  Xw  =  Total  stability  derivatives  of  the  longitudinal  X-force 
Z  =  Normal  force  along  the  body  Z-axis  (lbs) 

1,1  -  Total  stability  derivatives  of  the  normal  Z-force 

a  =  Remote  wind  angle  of  attack  relative  to  body  X-axis  (radians) 

ac  =  Rotor  angle  of  attack  (radians) 

Y  =  Lock  Intetia  number  =  — 

*b 

e  =  Qownwash  interference  angle  (radians) 
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V  sin  a  -  V. 

\  =  Rotor  inflow  ratio  =  - ^ - 

V  cos  a 

p  =  Rotor  tip  speed  ratio  =  - ^ - 

p  =  Air  density  (slug/cu  ft) 

a  -  Rotor  solidity  = 

0  =  Rotor  angular  velocity  (radians/second) 

0  =  Blade  twist  (radians) 

0  ^  =  Blade  section  pitch  angle  at  .75  radius  (radians) 

SUBSCRIPTS 

F  =  Pertaining  to  the  front  rotor 

FUS  =  Pertaining  to  the  fuselage 

T  =  Pertaining  to  the  horizontal  tailplane 

TR  =  Pertaining  to  the  tail  rotor. 

The  following  procedures  give  a  step-by-step  approach  that  can  be 
followed  to  find  the  desired  total  stability  derivatives.  The  calcu¬ 
lation  of  the  stability  derivatives  is  nothing  more  than  an  accounting 
of  the  forces  and  moments  acting  on  the  aircraft.  Because  of  the  number 
of  calculations  involved  and  the  interaction  of  many  different  compon¬ 
ents,  however,  the  calculation  procedure  is  necessarily  a  lengthy  one. 

C.  HELICOPTER  DATA 

Information  about  the  physical  layout  of  the  helicopter  must  be 
available  before  the  computations  can  begin.  The  following  must  be 
known: 
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a,  b,  c,  1,  R,  Q,  W,  1^,  Iy,  M^,  ^X’ 


(Cq  )  ,  (AR)y»  Sy,  3y 
0  T 

Because  of  the  high  aspect  ratio  of  a  helicopter  blade,  the  lift  curve 
slope  of  the  blade,  a,  can  generally  be  taken  as  5.73  per  radian.  Es¬ 
tablish  the  altitude  and  flight  speed,  V0,  at  which  the  derivatives  are 
desired  to  be  known.  The  altitude  will  determine  the  atmospheric  den¬ 
sity.  For  sea  level,  p  =  .002377  slugs/cu  ft. 

Knowing  the  above,  the  following  parameters  can  be  calculated: 

o,  p,  My,  V<-,  q0>  T.F.,  Rotor  Tip  Speed 


0.  CALCULATION  OF  HELICOPTER  TRIM  VALUES 

First  approximations  of  the  fuselage  lift  and  drag  coefficients  are 
determined  with  the  assumption  that  c*pUS  =  0.  C^  and  CQ  will  be  recal¬ 
culated  when  more  precise  information  is  known  about  <^<5.  A  theore¬ 
tical  or  experimentally  obtained  graph  of  C.  ,  Cn  and  CM  versus 

lFUS  FUS  FUS 

aFUS  1S  needed  t0  be91n  the  computation.  Figure  10.1-3  [Ref.  5]  is  an 
example  of  such  a  graph. 

The  first  approximation  of  fuselage  lift  and  drag  will  be  made  using 
the  coefficients  just  obtained. 


*-FUS  =  ^L  do  ^Z 
tud  i_F(JS  ^Fus 


(3-1) 


do 


(3-2) 
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The  initial  estimate  of  main  rotor  lift  and  drag  forces  can  now  be  made. 


L  =  W  - 


LFUS 


(3-3) 


0  = 


-D 


FUS 


(3-4) 


The  rotor  lift  and  drag  coefficients  are  non-dimensional  and  are 
obtained  by  dividing  main  rotor  lift  by  the  thrust  factor.  Both  are 
normalized  by  rotor  solidity. 


fe.)  =  _L_ 

\  a  /  T.F.o 


Reference  6  contains  theoretical  rotor  data  performance  charts  for 
hinged  rectangular  planform  blades  of  various  degrees  of  twist.  Differ¬ 
ent  charts  are  used  for  the  variables  of  blade  twist,  tip  speed  ratio 
and  advancing  tip  mach  number.  All  charts  are  based  on  a  rotor  solidity 
of  0.1.  Corrections  must  be  made  if  the  rotor  solidity  of  the  heli¬ 
copter  under  evaluation  varies  from  0.1.  The  charts  are  entered  knowing: 

r  1  r  1  r 

D  L  0 

— -  and  -  to  find  ,  a.  ,  a  and  6.75  . 

a  o  a  c 


Aa  =  a  -  0.1 


(3-7) 


(3-8) 


(3-9) 
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The  above  two  quantities,  \~^~j  g  i  and  \  o~J  i»  t09ether  with  the 
known  or  calculated  quantities  blade  twist  angle,  9,  main  rotor  tip 
speed  ratio,  p,  and  mach  number  of  the  advancing  blade  tip,  M^,  are 
needed  to  use  the  charts  in  [Ref.  5].  These  charts  will  yield  an  ap¬ 
proximation  of  the  following  main  rotor  trim  parameters: 


(1)  ac 

(2)  a0 

(3)  at 

(4)  bi 

(5)  9  75 

(6)  X 

(7) ^9 

a 


Rotor  angle  of  attack 

Blade  coning  angle 

Longitudinal  flapping  angle 

Lateral  flapping  angle 

Blade  section  pitch  angle  at  .75  radians 

Inflow  ratio 

Coefficient  of  rotor  torque  =  ^  ~p"r 


The  angle  of  attack  of  the  main  rotor,  ac,  and  the  rotor  torque,  Q, 
may  now  be  calculated  in  the  following  manner: 


a 


c 


(0fC)0. 


(3-10) 


Q  =  (T.F.)  aR 


(3-11) 


Interference  effects  can  exist  when  the  downwash  created  by  one 
aerodynamic  component  affects  the  performance  of  another  aerodynamic 
component.  These  effects  can  be  described  as  changes  in  the  local  angle 
of  attack  and  the  local  velocity.  Changes  in  the  local  velocity  are 
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usually  small,  however,  and  will  not  be  considered  here.  The  local 
angle  of  attack  can  be  expressed  by  the  following  equation: 

a LOCAL  =  “  +  1  ‘  £  (3'12) 

where  a  is  the  remote  wind  angle  of  attack  with  respect  to  the  X-axis 

i  is  the  geometric  inclination  of  the  aerodynamic  component  being 
evaluated  with  respect  to  the  X-axis 

e  is  the  aerodynamic  interference  angle. 

For  helicopters,  aerodynamic  interference  is  produced  mainly  by  the 
main  rotor  downwash.  The  downwash  velocity  of  a  rotor  varies  with  time 
as  well  as  position  and  the  result  is  an  exceedingly  complicated  situa¬ 
tion  to  evaluate.  Fortunately,  measurements  of  lift  and  drag  for  a 
single  rotor  helicopter  as  reported  in  Reference  5  show  that  inter¬ 
ference  effects  between  the  main  rotor  and  fuselage  and  the  main  rotor 
and  horizontal  tail  are  negligible. 

The  interference  between  the  front  and  rear  rotors  of  a  tandem  rotor 
configuration  is  more  significant.  Tandem  rotors  will  not  be  evaluated 
here,  however,  Reference  5  contains  charts  and  procedures  for  calcu¬ 
lating  this  interference  effect. 

The  downwash  interference  angles,  £,  are  usually  small.  They  can  be 
calculated  using  equation  (3-13). 

£FUS  =  eT  =  eTR  =  K  ('TAN  ac  ’  p  }  (3-13) 

where  K  is  the  rotor  interference  factor  and  is  usually  equal  to  1.0  for 
single  rotor  aircraft. 
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The  above  calculations  have  revealed  additional  information  about 

the  trim  condition  of  the  helicopter.  A  second  approximation  of  the 

fuselage  trim  angle  of  attack  can  now  be  made.  First  a  mathematical 

relationship  is  obtained  between  aPIK  and  CM  .  This  information  is 

rui  m  Fus 

then  plotted  against  experimental  data  for  the  same  two  factors.  The 

intersection  of  the  two  curves  will  yield  the  trim  value  for 

Use  equation  (3-14)  to  find  the  relation  between  ofp||C.  and  CM 

t-ub  mfus 


aFUS  ~ 


ebfi2M 

Dx  L  ‘  lz  D  +  - (at  +  ac  -  i  )]p  +  q0  Dx  S  a  (i  -  e)] 

-  D)  +  lz  L)p  -  q0  (lx  Sa)T 


'FUS 


do 


'FUS 


FUS 


-  <>X 


+  1. 


t)p  *  do  (lx  3a)j 


-  e 


FUS 


(3-14) 


The  only  two  unknowns  in  equation  (3-14)  are  tfpus  and  .  The 

suit  will  be  that  some  constant  times  Cu  will  equal  ar,,r .  Plot  this 

MPUs  FUS 


against  the  experimental  data  for  the  values  such  as  shown  in  Figure 
10.1-11  [Ref.  5].  The  intersection  of  the  two  lines  will  be  the  new 
approximation  of  cfp^^. 

The  recently  obtained  value  of  c<py<-  can  be  used  to  gain  a  second 

approximation  of  C,  ,  Cn  ,  CN  ,  and  CM  by  using  the  M,-||C|  vs. 

'“FUS  UFUS  FUS  FUS  ^ 

aFUS  ^1-9ure  a9ain- 

A  refined  approximation  of  the  fuselage  forces  can  be  gained  by 
using  equations  (3-1)  and  (3-2).  The  fuselage  pitching  moment  can  be 
calculated  from  equation  (3-15). 
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^FUS  =  AX  1 FUS 

rui  npys  aF(JS 


(3-15) 


Tail  Rotor  Calculations 

The  yawing  moment  of  the  fuselage  is  calculated  next. 


nfus  *  cnfus  \us  'fus 


(3-16) 


The  tail  rotor  thrust  must  be  equal  to  the  yaw  moments  produced  by 
the  main  rotor  and  the  fuselage.  Thus: 


T  _  nfus  qf 

1 TR  "  -  1Y 

XTR 


(3-17) 


The  coefficient  of  tail  rotor  thrust  can  now  be  obtained. 


— )  =  \—L 

o/TR  [(T.F. 


)CT 


(3-18) 


Knowing  the  values  of  0,  TR,  p.^,  and  Nl-j.  (calculated  previously)  plus 

C.  ’  TR 

and  assuming  (a  )-j.^  =  0,  use  the  appropriate  charts  in  [Ref.  5] 


TR 

to  find: 


(— )1  ■  (^) 

'  O  /  Tn  n  '  O  1 


T  ,  ,  ,\TR  ,  (0.  75)tr 

■J TR  o.l  '  J  TR 


compute: 


(V)  .  [(V)  .  *  (V) 

'  °  '  TR  l'  <J  /  q  1  2 p7  \  a  / 


(3-19) 


TR 
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The  drag  of  the  tail  rotor  can  now  be  calculated: 


>TR  =  Rr 


The  only  tail  rotor  derivative  needed  is: 


(3-20) 


3DTR  _  /T.F.  a\  iB  a  ) 

3uTD  \  0  R  /  Trj  \  dy  / 


(3-20a) 


The  last  term  in  equation  (3-20a)  is  obtained  from  the  appropriate 


/CL'\ 

chart  of  Reference  5  knowing  p,  My,  and  \-^-f 


The  characteristics  of  the  horizontal  tail  plane  are  determined  from 


the  following  equations: 


01  aFUS  +  eFUS 


Ofy  =  a  +  Ly  -  Cy 


CLy  ~  aj 


,(c  .cjL\ 

\  Do  "AR  ) 


Lr  ~  CLy  q°  ST 


°T  =  ^Dy  q°  ST 


(3-21) 


(3-22) 


(3-23) 


(3-24) 


(3-25) 


(3-26) 
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Forces  are  summed  in  the  X  and  Z  directions  to  obtain  a  new  approxi¬ 
mation  of  the  main  rotor  lift  and  drag.  The  following  equations  are 
used: 

Ki  =  W a  ’  t-pu$  (°f"eFUS^  ’  LT  ^a~eT^  +  °FUS  +  °T  *  °TR 

Kz  =  DFUS  ^a"eFUS^  +  °T  ^a_eT^  +  DTR  ^a"eTR^  +  LFUS  +  LT  ’  W 

where  Ki  represents  the  total  drag  being  developed  by  the  helicopter  and 
K2  represents  the  total  lift. 

The  lift  of  the  main  rotor  can  now  be  approximated  from  the  follow¬ 
ing  two  equations: 


(3-27) 

(3-28) 


Lc  =  Kx  a 


(3-29) 


La  -  K, 


°F  = 


(3-30) 


The  angle  of  attack  used  here  is  from  equation  (3-16). 

Equations  (3-5  through  3-9)  are  used  to  obtain  better  estimates  for 


(¥)■  (¥).,■  (¥)■“(¥> 


0.1 


The  steps  following  equation  (3-9)  to  equation  (3-25)  should  be 
repeated  until  the  trim  values  converge.  This  should  only  take  one  or 
two  iterations.  The  result  of  the  above  calculations  will  yield  the 
final  trim  values  of  the  helicopter  at  one  airspeed.  These  values  are 
needed  for  the  further  calculation  of  the  stability  derivatives. 
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E.  CALCULATION  OF  ROTOR  ISOLATED  DERIVATIVES 


Isolated  rotor  derivatives  are  defined  as  those  aerodynamic  para¬ 
meters  for  the  rotor  that  change  with  respect  to  tip  speed  ratio,  p, 
rotor  angle  of  attack,  ac,  and  blade  section  pitch  angle  at  0.75  radius, 
0  7g.  Theoretical  values  of  these  derivatives  are  plotted  on  charts  and 
presented  in  Reference  5  as  functions  of  cr,  9,  y,  and  p. 

These  charts  were  derived  for  rotor  solidity  equal  to  0.1.  If  the 
actual  rotor  solidity  for  the  helicopter  being  evaluated  differs  from 
0.1,  then  corrections  should  be  made  to  the  values  obtained  from  these 
charts. 

The  parameters  which  must  be  corrected  for  solidity  are  those  deal¬ 
ing  with  p  and  ac-  The  correction  factors  and  equations  to  use  when 
correcting  for  solidity  are  listed  below. 


Solidity  Corrections  for  (p)  Derivatives 


where 


Act  =  a  -  0 . 1 


(3-31) 


(3-32) 


1  +  Act 

C  1  " 
3(— ) 

CT 

da 

L  c  . 

(3-33) 


(  )q  .  -  denotes  stability  derivatives  for  rotor  solidity  a  =  0.1. 

These  values  can  be  directly  obtained  from  the  charts  of 
Section  7.5.  in  Reference  5. 
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Obtain  the  following  p  derivatives  from  Section  7.5  of  [Ref.  5]: 


3A 

3M 


3ai 


3p 


3M 


and  the  ac  derivatives: 


3ac  ’  3«c 


3  K  3a, 

da  '  3a 
c  c 


F.  CALCULATION  OF  THE  LOCAL  DERIVATIVES 

Using  the  values  obtained  from  the  charts  in  Reference  5  and  cor¬ 
rected  for  solidity,  if  necessary,  calculate  the  following  dimensional 
derivatives  for  the  main  rotor.  These  local  derivatives  are  dimensional 
expressions  of  the  change  of  local  forces  or  moments  of  various  compon- 
ents  with  respect  to  the  local  wind  conditions. 


u  derivatives: 


(T.F.)ct 

i-  L  -i 

9(-r) 

3Up 

a  R  J 

L  3p  J 

3Dp 

(7. F. )crl 

3u^ 

Q5  J 

L  9p  J 

3a1(: 

w 

3a  ,1 

3uF 

L  3mJ 

3Mhubf 

M  [ 

dalFl 

3uf 

L  2  j  1 

.aJJ 

lb-sec 

ft 


(3-42) 


lb-sec 

ft 


(3-43) 


rad-sec 

ft 


(3-44) 


lb-sec  (3-45) 
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ap  derivatives: 

3Lp 

[(T.F.)a] 

lb/rad 

3a  p 

.  a<*c  J 

30p 

C(T. F. )o]  | 

«¥)1 

Ib/rad 

3a  p 

L  *»c  J 

3MHUBf 

3a  p 

= 

[eb^] 

w 

lb- ft 
rad 

q  derivatives: 

3alp 

r  34 

i 

1 

07 

jO 

*n 

Lyoci. es3  - 

fjJ 

sec 

The  isolated  derivatives 

for 

the  fuselage  are  found 

next. 

figure  referred 

to  in  Section  D  and  the  trim  values 

for  the 

determine: 


3C 

3a 


lFUS 

FUS 


3Cr 


rUS 


3a 


FUS 


3C. 


'FUS 


3a 


FUS 


These  values  are  used  next  to  compute  the  local  derivatives  of 
lage.  These  are  dimensional  derivatives. 


UFUS  derivatives: 


9LFUS 

_  _2_ 

lfus 

lb-sec 

9ufus 

Vo 

ft 

a0FUS 

2 

dfus 

lb-sec 

3uFUS 

v0 

ft 

3Mfus 

_  _2_ 

M 

FUS 

lb-sec 

auFUS 

Vq 

ft 

(3-46) 

(3-47) 

(3-48) 

(3-49) 

Using  the 
fuselage, 

the  fuse- 

(3-50) 

(3*51) 

(3-52) 
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cfp^  derivatives: 

aLFUS  _ 
s“fus 

do 

*FUS 

(3Clfus) 
'3o,fus  ' 

1b 

rad 

(3-53) 

30fus  . 
9c,fus 

Qo  Ay 

afus 

(“"pus) 

V 9aFUS  ' 

1b 

rad 

(3-54) 

aMFUS  . 
3aFUS 

do  Ay 
*FUS 

(“Sus) 

'3aFUS  ’ 

1b- ft 
rad 

(3-55) 

The  derivatives 

for  the  horizontal 

tail  are  obtained 

using  the 

charts  and  procedures  in  Reference  7.  Determination  of  the 

local  de- 

rivatives  for  the 

horizontal 

tail  can 

be  made  using  the 

fol lowing 

equations: 

u  derivatives: 

3Ll 

3Uj 

=  h  lt 

lb-sec 

ft 

(3-56) 

30t 

-  2  n 

lb-sec 

3Uy 

'  V0  °T 

ft 

(3-57) 

a  derivatives: 

3Lt 

lb 

rad 

o?| 

— ?  1  - 
II 

do  Sy 

(3-58) 

1 

H- 

o 

(O 

2  lt 

lb 

(3-59) 

30fy 

tt( AR)y  aT 

rad 
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G.  CALCULATION  OF  THE  TOTAL  STABILITY  DERIVATIVES 


The  total  stability  derivatives  can  finally  be  determined.  The 
calculations  are  based  on  the  isolated  and  local  derivatives  found  in 
the  previous  sections. 

Some  angular  derivatives  needed  for  the  calculation  of  the  total 
stability  derivatives  are  as  follows: 


3aFUS 

3a 

3aj 

3a" 

1  "  Kft 

(*-+%) 

(3-60) 

3aFUS 

3u 

3aj 

3u" 

.  !et 

v0 

(  K  3\  \ 

\  M  3  H  / 

(3-61) 

3a 

3u 

krf(t 

3\\ 
3p  / 

R 

(3-62) 

3ap 

j 

1  "  krf 

(l-  — 
'  M 

3A  \ 

3aJ 
c  R 

(3-63) 

The  X- 

Force  Derivatives 

(a) 

X  (X 
u  u 

)  +  (X  )  ri  if  +  (X  ) 

P  v  u  FUS  u' 

T 

(3-64) 

where 

(X  ) 

U  F 

3LF 

3Up  a 

3Dp 

3uc 

(3-65) 
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„  3“fus 

3u 

3Lfus 

L3aFUS 
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1 


(b) 


(3-68) 


I 


(b) 
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The  Pitching  Moment  (M)  Derivatives 
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where 


(3-84) 


(3-85) 


(3-86) 
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IV.  SOLUTION  OF  THE  CHARACTERISTIC  EQUATION 

In  Chapter  II  it  was  shown  that  the  longitudinal  dynamics  of  the 
helicopter's  motion  could  be  described  by  three  equations.  These  equa¬ 
tions  deal  with  forces  in  the  X  and  Z  directions  and  a  moment,  M,  about 
the  pitch  axis.  Knowing  the  stability  derivatives,  as  outlined  in 
Chapter  III,  the  stage  is  set  for  determining  the  motions  that  charac¬ 
terize  the  helicopter's  response  to  these  forces  and  moments. 

The  pro  2dure  for  obtaining  the  aircraft's  modes  of  motion  lies  in  a 
simultaneous  solution  to  the  three  equations  for  longitudinal  motion.  A 
determinant  is  ideally  suited  for  this  purpose.  What  is  sought  is  the 
aircraft's  response  to  gravitational  and  aerodynamic  forces  and  moments. 
Therefore,  the  stability  derivatives  are  evaluated  and  the  three  equa¬ 
tions  set  equal  to  zero  and  solved. 

The  determinant  in  question  is  derived  from  the  two  force  equations 
and  one  moment  equation  that  have  been  used  for  longitudinal  dynamics. 
For  the  case  of  studying  the  helicopter's  natural  response  (i.e.,  no 
forcing  functions)  the  determinant  is  set  equal  to  zero  and  the  modes  of 
motion  are  then  found.  If  the  helicopter's  response  to  a  given  input  is 
desired,  then  the  determinant  is  set  equal  to  that  input  (the  forcing 
function)  and  solved.  An  example  of  such  a  forcing  function  is  the  step 
input  of  cyclic  or  collective  movement. 

The  determinant  derivation,  from  Reference  6,  starts  with  the  lin¬ 


earized  equations  of  motion: 


^  u  =  -  w0  COS  X  +  AX  (4-1) 

^w-^jv0  =  -w0sinx+AZ  (4-2) 

p  0  =  AM  (4-3) 

AX  and  AZ  are  the  forces  arising  from  aerodynamics  in  distrubed 
flight  and  M  is  the  aerodynamic  moment  in  disturbed  flight.  0  =  q. 

It  was  shown  in  Chapter  II  that  the  changes  in  forces  and  moments 
arising  from  small  disturbances  could  be  written  as  Taylor  series  ex¬ 
pansions.  Thus  a  substitution  for  AX  can  be  made  as, 


_  ax  sx  ^  ax  ^  ax  D  _  ax  „ 

iX  =  53  u  *  55  u  *  55  *  5bT  Bl  *  3^  e» 


(4-4) 


where  3X  and  0O  are  cyclic  and  collective  pitch  control  terms, 
respecti vely. 

3  X 

Using  a  shorthand  notation  of  |-r  =  Xi ,  equation  (4-4)  can  be 
written  as: 

AX  =  X  u  +  X  w  +  X  q  +  XD  B,  +  X_  0O  (4-5) 

u  w  q  ^  Bx  1  0O 


The  linearized  equations  (4-1),  (4-2),  and  (4-3)  can  now  be  written 
as: 

-  u  -  X  u  -  X  w-X  q  +  w0  cos  x  =  Xn  Bt  +  X  0O  (4-6) 

g  u  w  q  M  Bi  1  0O  0 

'  Zu  U  +  5  "  _  Zw  w  '  Zq  q  ‘  i  V®  +  w6  Sin  1  =  ZS1  +  Z0O  00  (4"7) 

-  M  u  -  M  w-M*w+B0-M  q  =  MD  Bi  +  M  0q  (4-8) 

u  w  w  q  M  Bi  1  90 
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The  resulting  equations  are  linear  ones  with  constant  coefficients, 
and  can  be  written  in  determinant  form  with  Laplacian  notation  for 
compactness.  With  this  notation  S  represents  d(  )/dt.  The  resulting 
determinant  is  set  equal  to  zero  when  the  natural  response  of  the  air¬ 
craft  is  to  be  examined.  (Bi  and  0O,  the  cyclic  and  collective  pitch 
inputs  =  zero.)  It  should  be  noted  that  2^  is  always  equal  to  zero  and 
that  and  Xq  are  negligibly  small  so  that  the  final  form  of  the 
stability  determinant  is: 


(4-9) 


Note:  z  is  usually  a  small  angle  such  that  cos  t  =  1.  For  com¬ 

patibility  of  units,  velocity  must  be  in  feet  per  second.  To  convert 
knots  to  FPS,  multiply  by  1.6889. 
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A.  THE  CHARACTERISTIC  EQUATION 

The  results  of  the  determinant  generated  in  the  above  manner  is 
called  the  characteristic  equation  and  in  general  is  a  quartic  equation. 
It  will  be  shown  that  in  a  hover  the  characteristic  equation  is  reduced 
to  a  cubic  equation. 

The  characteristic  equation  has  the  following  form: 


AS4  +  BS3  +  CS2  +  DS  +  E  =  0 


(4-10) 


where  the  value  of  the  coefficients,  A  thru  E,  are  determined  by  the 
determinant  cross-product. 
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Solving  the  characteristic  equation  will  yield  four  roots.  In 
classical  fixed-wing  aircraft  longitudinal  analysis,  the  solution  ac¬ 
tually  yields  two  sets  of  complex  roots.  These  sets  of  roots  describe 
two  sinusoidal  motions  of  the  aircraft  which  are  distinguished  from  each 
other  by  their  periods.  The  longer  period  motion  is  called  the  phugoid 
and  the  shorter  period  motion  is  called  the  short  period. 

The  fixed-wing  characteristic  equation  quartic  can  also  be  repre¬ 
sented  by  the  product  of  two  quadratic  equations.  The  coefficients  of 
the  two  quadratic  equations  contain  terms  which  define  the  damping 
ratio,  £,  and  the  response  frequency,  u>,  of  the  modes  of  motion  associ¬ 
ated  with  the  quadratic  equations. 

AS4  +  BS3  +  CS2  +  DS  +  E  =  js2  +  2CpwpS  +  Wp2  j 

Is2  +  ^SP'V5  +  msp|  '  (4"n) 

where  the  subscript  P  defines  that  motion  associated  with  the  long 
period  motion  of  the  aircraft,  or  phugoid,  and  the  subscript  SP  it  used 
to  denote  the  short  period  motion. 

The  phugoid  mode  of  motion  is  one  in  which  the  aircraft's  angle  of 
attack  remains  essentially  constant  while  airspeed  and  altitude  change 
as  aircraft  kinetic  energy  (airspeed)  and  potential  energy  (altitude) 
are  exchanged  until  the  aircraft's  motion  dampens  out  at  trim  airspeed 
(when  the  system  is  convergent)  or  until  the  aircraft  departs  controlled 
flight  (if  the  phugoid  has  a  divergent  nature)  (see  Figure  4-1).  The 
respective  velocities  and  altitudes  at  points  1,  2,  and  3  would  be  in¬ 
creasing  or  decreasing  depending  on  whether  the  aircraft's  oscillation 
was  damped  or  divergent  as  shown  in  Figure  1-2. 
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Figure  4-1.  Phugoid  Response 


For  helicopters  the  period  of  the  phugoid  is  typically  long,  on  the 
order  of  30  seconds  while  the  short  period  mode  of  motion  is  a  heavily 
damped  one  in  which  the  aircraft's  velocity  remains  constant.  Since  the 
short  period  is  heavily  damped,  its  period  is  typically  on  the  order  of 
one  second. 

The  aircraft's  resultant  motion  can  be  compared  to  a  spring-mass- 
damper  system  which  has  two  springs  and  two  dampers.  Keeping  this  in 
mind,  £,  the  damping  ratio,  is  a  measure  of  the  amount  of  damping  in  the 
system.  Insofar  as  the  short  period  has  a  more  heavily  damped  motion, 
it  follows  that  the  value  of  Csp  will  be  higher  than  that  for  Cp. 
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B.  HELICOPTER  RESPONSE  IN  A  HOVER 


Unfortunately,  the  helicopter's  modal  response  cannot  be  given  the 
same  clear,  consistent  physical  interpretation  as  that  for  a  fixed-wing 
aircraft.  One  reason  for  this  becomes  immediately  apparent  when  hover¬ 
ing  flight  is  considered.  For  a  conventional  fixed-wing  aircraft  all 
phases  of  flight  are  conducted  with  some  forward  flight  velocity.  This 
is  obviously  not  so  with  the  helicopter. 

Suprisingly,  the  characteristic  equation  for  a  helicopter  in  hover 
is  a  cubic  equation,  not  a  quartic.  The  reason  for  this  is  found  by 
examining  the  Z-force  equation  in  the  determinant: 

Zu  S-Zw  V0S  =  0  (4-12) 

V0  equals  zero  in  a  hover,  and  if  Z^  equals  zero  for  hovering  flight 
(as  it  frequently  does)  the  vertical  motion  is  immediately  decoupled 
from  the  pitching  motions  and  fore  and  aft  motions.  As  a  result,  the 
vertical  motion  is  entirely  dependent  on  Z and  is  usually  a  heavily 
damped  subsidance  motion  since  Z^  is  usually  a  large  negative  number. 
The  assertion  that  the  vertical  motion  of  the  helicopter  in  a  hover  is 
very  dependent  on  Z^  is  borne  out  by  the  fact  that  a  helicopter's 
vertical  motion  is  known  to  be  very  responsive  to  vertical  gusts  of  wind 
and  to  collective  pitch  inputs. 

Once  the  solution  to  the  Z-force  equation  is  known  (S  =  Z  ) ,  it  can 
be  removed  from  the  determinant  and  the  system  of  equations  can  then  be 
reduced  to  a  2  'A  2  determinant. 
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1 


S-Xu  w 

-M  S2-M  S 

u  q 


=  0 


and  the  resulting  characteristic  equation  will  be  a  cubic: 


(4-13) 


S3  - 


(X, 


+  Mq>  S2  * 


X  M  S  + 
u  q 


MuW  =  0 


(4-14) 


The  solution  of  the  cubic  equation  will  generally  yield  a  negative 
real  root  and  a  positive  complex  root  pair.  The  negative  real  root 
indicates  a  stable  convergence  and  is  principally  due  to  pitch  damping 
of  the  rotor  whereas  the  positive  pair  of  complex  roots  indicates  un¬ 
stable  oscillation  and  is  due  to  the  coupling  of  the  pitch  and  longi¬ 
tudinal  velocity  by  the  speed  stability  derivative  M^. 

For  an  articulated  rotor,  the  real  root  (which  indicates  the  short  < 

period  mode)  typically  has  a  time  to  half  amplitude  of  Th  =  1  to  2  i 

| 

seconds.  The  long  period  mode  is  represented  by  the  oscillatory  root  ! 

i 

I 

and  has  a  period  of  T  =  10  to  20  seconds  and  since  it  is  a  divergent 
motion,  a  time  to  double  amplitude  of  T2  =  3  to  4  seconds  (Ref.  10). 

Although  the  phugoid  motion  is  unstable,  the  period  and  time  to 
double  amplitude  are  sufficiently  long  for  the  pilot  to  observe  the 
helicopter's  reactions  and  make  necessary  control  movements  to  maintain 
control  of  the  aircraft. 

Hingeless  rotors  have  a  higher  degree  of  pitch  damping  than  articu¬ 
lated  rotors.  This  high  degree  of  damping  serves  to  greatly  increase 
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the  magnitude  of  the  real  root  and  it  also  increases  the  period  and  time 
to  double  amplitude  of  the  oscillatory  mode.  For  hingeless  rotors  the 
time  to  half  amplitude  is  typically  0.2  to  0.5  seconds,  while  the 
oscillatory  phugoid  mode  has  a  period  of  10  to  20  seconds  with  a  time  to 
double  amplitude  of  10  to  15  seconds. 


C.  APPROXIMATION  OF  THE  SHORT  PERIOD  DURING  HOVER 

The  initial  response  of  the  helicopter  to  gusts  is  primarily  that  of 
vertical  and  pitch  acceleration,  with  little  longitudinal  acceleration. 
Since  this  analysis  assumes  perturbations  from  steady,  level  flight,  the 
longitudinal  degree  of  freedom  can  be  neglected  in  order  to  approximate 
the  short  period  mode.  (It  should  be  noted  that  the  short  period  mode 
is  heavily  damped  and  is  characterized  by  near  zero  velocity  change,  so 
that  the  neglect  of  longitudinal  acceleration  is  reasonable  for  this 
analysis. ) 

With  this  assumption,  the  stability  determinant  is  reduced  to: 


-M 

w 

s-z. 


w 


S-M. 


-V 


=  0 


and  the  characteristic  equation  becomes: 


(4-15) 


(S  -  M  )  (S-Z  )  -  M  V  =  0  (4-16) 

q  w  w  v 

It  has  been  explained  that  in  hover  the  pitch  and  vertical  motions 
decouple  and  the  two  solutions  of  the  characteristic  equations  are: 


S  =  Z  and  S  =  M 

w  q 


(4-17) 
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As  a  check  on  the  validity  of  this  assertion  note  the  following 
comparison  in  Table  4-1  between  stability  derivatives  and  roots  from  the 
solution  of  the  stability  quartic  for  three  different  helicopters.  The 
comparison  shows  generally  good  correlation. 


TABLE  4-1.  Comparison 


— 

Zw 

hq 

ROOTS 

BO- 105 

-0.3317 

-3.3972 

-0.331,  -3.4521 

CH-53D 

-0.298 

-0.499 

-0.2934,  -0.8232 

OH-6 

-0.3404 

-1.7645 

-0.3544,  -1.8794 

The  helicopter  has  neutral  static  stability  when  disturbed  by  per¬ 
turbations  in  pitch  or  roll.  This  is  because  no  moments  are  generated 
directly  by  these  motions  to  move  the  helicopter  back  to  or  away  from 
its  equilibrium  position.  However,  the  helicopter  in  hover  does  possess 
positive  static  stability  when  disturbed  by  longitudinal  or  lateral 
perturbations  of  wind  velocity. 

D.  HELICOPTER  RESPONSE  IN  FORWARD  FLIGHT 

As  the  helicopter  departs  hovering  flight  and  transitions  to  forward 
flight,  new  forces  come  into  existance  which  start  to  change  the  stabil¬ 
ity  picture.  These  forces  arise  from  the  increasing  dynamic  pressure 
being  built  up  on  the  aerodynamic  surfaces  of  the  helicopter.  As  a 
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result  the  dynamics  of  the  helicopter  in  forward  flight  are  different 
from  those  characteristics  exhibited  in  hover.  The  forces  and  moments 
that  act  on  the  helicopter  are  contributed  by  the  main  rotor,  tail 
rotor,  fuselage  and  tail  surface  aerodynamics,  and  gravity. 

A  stability  derivative  that  arises  from  forward  flight  is  M  ,  the 
pitching  moment  due  to  angle  of  attack  perturbation.  This  derivative  is 
usually  equal  to  zero  for  the  hover  condition.  It  has  an  unstabilizing 
effect  on  the  helicopter's  motion  for  the  following  reason:  As  the 
helicopter  accelerates,  the  advancing  blade  experiences  increased 
dynamic  pressure.  This  increase  in  pressure  results  in  a  change  in  the 
blade  angle  of  attack  which  produces  a  lateral  moment  on  the  rotor  disc 
(toward  the  retreating  blade).  The  moment  thus  generated  is  propor¬ 
tional  to  forward  velocity.  Because  of  the  gyroscopic  effect  of  the 
rotating  disc,  the  tip  path  plane  will  respond  to  this  moment  after  90 
degrees  of  movement.  Therefore  the  rotor  disc  will  be  tilted  aft.  The 
angle  of  attack  increase  results  in  a  pitch  up  moment  of  the  aircraft 
which  further  increases  the  angle  of  attack  on  the  rotor  system.  There¬ 
fore  the  dynamics  of  the  rotor  are  a  source  of  instability  for  the 
helicopter  in  forward  flight. 

To  counter  the  angle  of  attack  instability  of  the  main  rotor  in  for¬ 
ward  flight,  a  horizontal  tail  can  be  incorporated  on  the  helicopter. 
The  forces  and  moments  produced  by  the  horizontal  tail  are  proportional 
to  forward  velocity  so  that  they  are  approximately  zero  in  hover  and 
increase  with  speed.  The  horizontal  tail  will  have  a  stabilizing  in¬ 
fluence  on  the  helicopter's  motions  in  much  the  same  way  as  a  horizontal 
tail  on  a  fixed-wing  aircraft.  Since  the  moments  produced  by  both  the 


main  rotor  and  horizontal  tail  are  proportional  to  velocity,  their 
relative  contributions  to  stability  are  actually  independent  of  speed. 

As  can  be  seen,  Mw  is  influenced  by  two  opposing  sources,  an  un¬ 
stable  contribution  from  the  rotor  and  a  stabilizing  contribution  from 
the  horizontal  tail.  Without  a  horizontal  tail,  the  dynamics  of  the 
helicopter  in  forward  flight  are  characterized  by  two  stable  damped 
motions  (from  the  negative  real  roots)  and  an  unstable  oscillatory  mode 
(from  a  positive  pair  of  complex  roots).  In  this  configuration  the 
flying  qualities  are  degraded  due  to  the  angle  of  attack  instability. 

The  dynamics  of  the  helicopter  in  forward  flight  can  be  changed  by 
the  addition  of  a  large  enough  horizontal  tail  such  that  static 
stability  is  achieved.  For  the  case  of  a  large  horizontal  tail,  the 
pitch  and  vertical  real  roots  are  transformed  into  two  oscillatory  roots 
with  a  short  period  and  high  damping.  The  latter  motion  is  similar  to 
that  found  in  a  fixed-wing  aircraft. 

In  actual  practice  other  considerations  must  be  taken  into  account 
which  may  eclipse  the  goal  of  achieving  fixed-wing- 1  ike  dynamics  in 
forward  flight.  One  such  limitation  is  that  a  horizontal  tail  which  is 
large  enough  to  counter  the  unstabilizing  influences  of  the  main  rotor 
may  simply  be  too  large  for  weight  or  drag  considerations.  Another 
factor  to  consider  is  that  tail  effectiveness  is  reduced  at  low  speeds 
due  to  interference  with  rotor  and  fuselage  wakes.  In  spite  of  these 
problems  it  should  be  noted  that  almost  all  single  rotor  helicopters  do 
have  horizontal  tail  planes  because  of  the  improvement  in  flying  quali¬ 
ties  that  the  addition  of  this  component  provides. 


E.  APPROXIMATION  OF  THE  SHORT  PERIOD  IN  FORWARD  FLIGHT 

The  assumptions  made  for  the  analysis  of  the  short  period  motion  in 
a  hover  are  also  valid  for  the  approximation  of  the  short  period  motion 
in  forward  flight.  Of  course  the  characteristic  equation  (4-15)  is  not 
solved  so  readily  because  of  the  -MwV  term.  In  a  hover  this  term  went 
to  zero,  but  in  forward  flight  the  characteristic  equation  must  be 
solved. 


(S  *  Mq)  (S  -  Zw)  =  MwV  (4-18) 

Typically,  the  short  period  motion  of  the  helicopter  is  character¬ 
ized  by  two  negative  real  roots.  This,  of  course,  means  that  the  short 
period  is  a  stable  motion  and  is  heavily  damped.  Because  of  the  action 
of  the  horizontal  tail  as  explained  above,  certain  helicopters  do 
occasionally  exhibit  a  stable  oscillatory  mode,  more  in  keeping  with 
fixed-wing  dynamics. 
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INTERPRETATION  OF  THE  STABILITY  DERIVATIVES 


It  has  been  shown  how  the  individual  stability  derivatives  can  be 
calculated  knowing  some  basic  data  about  the  helicopter,  and  having 
access  to  theoretical  charts  relating  aircraft  performance  and  response 
parameters.  An  explanation  of  the  characteristics  of  these  derivatives 
and  their  effects  on  the  helicopter's  motion  is  now  in  order.  To  begin 
with  it  should  be  noted  that  since  all  airframe  contributions  are  pro¬ 
portional  to  airspeed,  the  fuselage  and  horizontal  tail  forces  will  be 
equal  to  zero  at  zero  flight  velocity.  This  fact  serves  to  simplify 
certain  calculations  for  the  hover  condition. 

A.  Xu  =  DRAG  DAMPING 

This  stability  derivative  acts  as  a  damping  force.  It  will  be 
negative  in  sign.  The  interpretation  of  X^  is  that  it  represents  an 
increase  in  drag  with  an  increase  in  forward  flight  velocity.  Physi¬ 
cally  this  is  seen  to  be  true  because  as  forward  speed  increases  the 
thrust  vector  (and  the  rotor  disk)  must  be  tilted  more  forward  to  over¬ 
come  the  effects  of  increased  drag. 

In  relation  to  the  dynamic  motions,  Xu  has  a  weak  effect  on  the 
phugoid,  but  one  which  tends  to  make  the  phugoid  more  stable. 
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B.  X  =  DRAG  DUE  TO  ANGLE  OF  ATTACK 

w 

Xw  can  be  interpreted  as  the  change  in  drag  on  the  aircraft  which  is 
brought  about  by  change  in  angle  of  attack.  It  is  usually  of  small 
value  and  does  not  have  much  influence  on  either  the  static  or  dynamic 
characteristics  of  stability. 

X^  is  usually  very  small  or  zero  for  hovering  flight.  For  the 
purposes  of  simplifying  calculations,  it  can  safely  be  assumed  to  equal 
zero  for  the  hover  condition. 

C.  Zu  =  LIFT  DUE  TO  VELOCITY 

This  stability  derivative  is  always  negative  for  fixed-wing  aircraft 

and  corresponds  to  increased  lift  at  higher  velocities.  (Remember,  the 

Z-axis  is  positive  downward.)  This  is  not  so  for  helicopters,  however. 

According  to  Reference  9,  for  helicopters  Zu  is  negative  at  low  speeds 

but  positive  at  high  speeds.  This  is  not  especially  significant  except 

at  higher  forward  velocity  where  it  might  affect  the  dynamic  divergence. 

Z  is  usually  small  for  an  articulated  rotor. 

Like  X  ,  Z  can  be  assumed  to  be  equal  to  zero  for  the  hover  con- 
w’  u  M 

dition  and  for  the  same  reasons. 

0.  Z  =  VERTICAL  QAMPING 
w 

Zw  acts  as  a  damping  force  in  the  same  manner  that  does.  It  also 
is  negative  in  sign  and  occurs  because  of  the  vertical  motions  of  the 
aircraft.  For  a  helicopter  this  is  an  important  parameter,  especially 
in  a  hover  where  it  describes  the  response  of  the  aircraft  to  vertical 
gusts.  This  stability  derivative  is  nearly  independent  of  airspeed  for 
the  helicopter  whereas  it  is  proportional  to  airspeed  for  fixed-wing 


aircraft. 


E.  Mu  =  VELOCITY  STABILITY 

This  parameter  describes  the  pitch  tendencies  of  the  aircraft  with 
respect  to  speed  changes.  According  to  Reference  9,  at  hover  and  for 
very  low  speeds,  most  helicopter  configurations  have  Mu  positive.  This 
leads  to  a  positive  stick  gradient.  A  positive  value  for  could  lead 
to  oscillatory  instability  and  it  indicates  the  aircraft  is  sensitive  to 
turbulence. 

changes  sign  at  high  forward  flight  velocities.  Negative  values 
of  lead  to  dynamic  divergence. 

F.  M  =  ANGLE  OF  ATTACK  STABILITY 

w 

indicates  the  tendency  of  the  aircraft  to  pitch  up  or  down  as 
angle  of  attack  is  increased.  A  negative  value  of  is  stabilizing  as 
it  tends  to  return  the  aircraft  to  its  previous  position  whereas  a 
positive  value  would  be  divergent  in  nature.  Most  helicopters  exhibit 
values  of  M  that  are  neutral  or  positive.  A  positive  value  will  lead 
to  dynamic  divergence  in  forward  flight.  Seckel  claims  that  the  center 
of  gravity  position  can  have  an  effect  on  the  values  of  and  conse¬ 
quently  on  the  resultant  stability  characteristics  of  the  aircraft.  If 
the  fuselage  and  horizontal  tail  contributions  are  stable  with  respect 
to  angle  of  attack,  moving  the  center  of  gravity  forward  will  make  the 
aircraft  more  stable  and  aft  less  so.  Conversely,  if  the  fuselage  and 
horizontal  tail  have  destabilizing  tenencies,  then  the  opposite  will 
occur. 

According  to  Reference  10,  -eceives  an  unstable  contribution  from 
the  rotor  and  fuselage  and  a  stable  contribution  from  the  horizontal 


tail.  M  is  the  third  stability  derivative  which  can  be  assumed  to  be 
w 

equal  to  zero  for  hover  calculations  with  no  loss  of  accuracy  in  the 
results. 

G.  Mq  =  PITCH  DAMPING 

This  derivative  is  negative  and  is  a  very  important  one  for  response 
to  control  deflection  and  dynamic  stability.  Most  helicopters  require 
some  augmentation  of  this  angular  damping  for  good  handling  qualities. 
The  numerical  value  of  Mq  should  be  less  than  -0.5  [Ref.  9].  Values  of 
-2.0  are  better,  but  the  degree  of  improvement  decreases  for  values  less 
than  -2.0. 
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VI.  ANALYSIS  OF  RESULTS 


As  can  be  seen  from  the  sample  problem,  Appendix  A,  the  stability 
characteristics  of  the  helicopter  under  investigation  consisted  of  a 
pair  of  damped  oscillatory  roots,  a  divergent  real  root,  and  a  conver¬ 
gent  real  root.  These  characteristics  differ  from  the  low-speed 
response  for  the  same  helicopter.  At  advance  ratios  of  0.1  to  0.2  this 
aircraft's  response  was  identical  to  that  expected  for  helicopters  in 
low  speed  flight,  namely  two  heavily  damped  stable  roots  and  a  pair  of 
divergent  oscillatory  roots.  Thus,  the  stability  characteristics  of  the 
helicopter  have  changed  at  high  speed. 

The  resulting  motion  of  the  helicopter  at  high  advance  ratios  cannot 
be  accurately  predicted,  however,  as  can  be  seen  from  Table  6-1.  Tnis 
table  shows  tne  different  response  modes  of  five  helicopters  at  various 
conditions  of  airspeed  and  center  of  gravity  location.  Some  interesting 
observations  can  be  made  from  this  table. 

The  response  of  some  helicopters  may  not  change  at  all  with  airspeed. 
The  CH-53D  is  an  example  of  this.  Data  obtained  from  Reference  11 
indicates  the  CH-53  retains  the  response  mode  characteristics  associated 
with  hovering  flight  out  to  at  least  140  knots.  Different  responses  are 
exhibited  by  other  helicopters.  The  OH-6  has  a  typical  hover  response 
mode  which  changes  into  two  damped  real  roots  and  3  damped  oscillatory 
motion  prior  to  100  knots.  The  30-105  also  changes  from  the  hover 
response  at  forward  speed  but  not  until  an  airspeed  greater  than  100 
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TABLE  6-1.  RESPONSE  MODE  COMPARISON 


AIRCRAFT 

ROOTS 

HOVER 

100  KNOTS 

140  KNOTS 

CH-53D 

0.1139  ±  0.4705 

0.1648  ±  0.3097 

0.2538  ±  0.2769 

-0.2934 

-0.3374 

-0.2793 

-0.8232 

.... 

-1.3822 

-1.7994 

HOVER 

100  KNOTS 

130  KNOTS 

OH-6 

0.0516  ±  0.4658 

-0.0192  ±  0.3186 

-0.1027  ±  0.3729 

-0.3544 

-1.0092 

-1.0033 

-1.8794 

-2.6212 

-2.8293 

HOVER 

100  KNOTS 

145  KNOTS 

BO- 105 

0.0188  ±  0.4333 

0.2709  z  0.4673 

1.2449 

-0.331 

-0.4040 

0.3999 

-3.4521 

-4.7375 

-0.3189 

-5.9681 

HOVER 

100  KNOTS 

130  KNOTS 

UH-1H 

(MID  C.G. ) 

0.1685  ±  0.3535 

-0.0060  ±  0.2S20 

-0.0211  ±  0.2392 

-0.4577  ±  0. 1527 

-0.3653  ±  1.0150 

-0.9354  ±  1.4882 

(AFT  C.G.  ) 

0.0220  t  0.2380 

-0. 3720  ±  0. 7592 

HCVER 

ICO  KNOTS 

140  KNOTS 

AH-  IG 

(AFT  C.G. ) 

0. 1205  t  0. 2517 

0.  0089  ±  0.2666 

-0.0337  ±  0.2051 

-0.4323  ±  3.1989 

-0.7179  t  0.5312 

-0.  '442  :  1.5398 
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knots  and  then  the  response  is  represented  by  four  real  roots,  two 
convergent  and  two  divergent. 

Still  other  examples  of  stability  can  be  seen  with  the  UH-1H  and  the 
AH-1G.  Both  of  these  helicopters  have  stable  oscillatory  periods 
typical  of  the  response  of  a  fixed-wing  airplane.  Of  interest  here  is 
the  change  in  response  with  change  in  center  of  gravity  position  at  100 
knots.  Moving  the  center  of  gravity  can  have  a  profound  effect  of  a 
helicopter.  For  both  the  UH-1H  and  the  AH-1G,  moving  the  center  of 
gravity  aft  has  a  destabilizing  influence. 
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APPENDIX  A 
SAMPLE  PROBLEM 


MAIN  ROTOR 

R  =  24  ft 
Q  =  29  rad/sec 
b  =  4 
e  =  0.5  ft 
0!  =  -8° 


FUSELAGE 

W 

:y 

1 FUS 


=  10,000  lbs 
=  17,500  slug- ft2 
=  50  ft 


TAIL  ROTOR 

R  =  4.6  ft 

b  =  4 

c  =0.75  ft 
0!  =  0° 


c  =1.75  ft 
a  =  5.73/rad 
M$  =85.4  slug- ft 
Ib  =  1200  slug- ft2 
17  =  -6  ft 

lv  =  -10  ft 
XF 

Ax  =  60  ft2 
Az  =  200  ft2 


aTR  =  5.73/rad 

1 v  =  -30  ft 
i  R 

17  =  -6  ft 

^TR 


i  =  0 


Q  =  146.6  rad/ sec 


HORIZONTAL  TAIL 

S  =  20  ft2  aT  =  3.5/rad 

AR  =  4.0  lv  =  -20  ft 

XT 


=  0 


=  0 


Airfoil  =  NACA  0015 
Cqo  =  0.02 


The  helicopter's  stability  characteristics  will  be  evaluated  at 
V0  =  203  FPS  and  at  sea  level.  Using  the  data  supplied  for  the  heli¬ 
copter,  and  the  definitions  listed  in  Chapter  III  or  the  program  STAB 
from  Appendix  8,  the  following  values  are  determined: 


=  0.002378  slugs/ft3 

T.F.  =  2,084,504  lbs 
m  =  310.8  slugs 
QR  =  696  fps 
V-j.jp  =  899  fps 

(T.F.)tr  =  71,896 

(S1R)tr  =  674  fps 

(M,)  =  0.783 

'  TR 


q0  =  49.0  psi 

o  =  0.093 

p  =  0.292 

Myip  =  0.805 

aTR  =  0.207 
Mtr  =  °-30 


Use  Figure  A-l  to  obtain  values  for  C,  and  Cn  for 

lFUS  uFUS 

C.  =  -0.005 
LFUS 


c 


D 


FUS 


0. 16 


0. 
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Equations  (3-1)  and  (3-2)  yield  the  fuselage  lift  and  drag: 


lfus  =  “49  lb 


Dfus  =  470.4  lb 


Main  rotor  lift  and  drag  are  obtained  from  equations  (3-3)  and  (3-4). 


LF  =  10049  lb 


Dp  =  -470.4  lb 


The  main  rotor  lift  and  drag  coefficients  are  calculated  from  equations 
(3-5)  and  (3-6). 


L  D 

—  =  0.0473  —  =  -0.0022 

a  a 

Equations  (3-7),  (3-8),  and  (3-9)  are  used  to  correct  these  values  if 
the  main  rotor  solidity  differs  from  0.1. 


“  =  -°- 007  ( V)0  J  =  °- 0473  (^-)0  J  -  -0.  0021 

Now  use  the  appropriate  charts  in  Reference  6  to  obtain  values  for  the 
fol lowing: 


at  =  0.0611  rad  ac  =  -0.1222  rad 

-2  =  0.0025  0  7,  =  5° 

o  .  75 

Also  use  the  charts  in  Section  5.3  of  Reference  5  to  obtain: 

\  =  -0.045  a0  =  2.3° 
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The  angle  of  attack  of  the  main  rotor  and  the  main  rotor  torque  are 


calculated  from  equations  (3-10)  and  (3-11). 

ac  =  -0.1241  rad  Q  =  11632  ft- lb 


The  downwash  interference  factors  for  this  configuration  are  Kppy,.  = 
Kp-j-  =  KpTR  =  1.0  and  the  downwash  interference  angles  equal  0.0312  rad¬ 
ians  from  equation  (3-13). 

A  relation  between  acllc  and  Cu  is  needed  to  plot  against  the 

FUS  Mfus 

experimentally  obtained  fuselage  pitching  moment  data.  Equation  (3-14) 
will  yield: 


ofruc  =  -0.0164  +  0.5911  CM 
FUS  "pus 

The  point  of  intersection  on  Figure  A-2  will  yield  the  fuselage  trim 
angle  of  attack,  -2  degrees  or  -0.0349  radians.  Using  this  value  use 
Figure  A-l  again  to  obtain  the  following  new  values  for  the  parameters 
i ndicated: 


C.  =  -0.0065 
LFUS 


-0.025 


Cn  =  0.158  Cw  =  0.007 

UFUS  FUS 

Use  equations  (3-1)  and  (3-2)  to  recalculate  Lpy<.  and  Dp^.  Also  use 
equations  (3-15)  and  (3-16)  to  determine  MpUS  and  Npy$,  respectively. 

LF(js  =  -63.7  lb  Mfus  =  -3675  ft- lb 

DF(js  =  464.5  lb  NpyS  =  1029  ft- lb 
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Using  Np^  and  Qp  determine  the  tail  rotor  thrust  and  the  tail  rotor 
lift  coefficient  from  equations  (3-17)  and  (3-18). 


TTR  =  422  lb 


0.0293 


Enter  the  charts  in  Reference  6  to  obtain  additional  tail  rotor  para¬ 
meters,  knowing  the  tail  rotor  lift  coefficient.  Blade  twist  =  0, 
p  =  0.3,  and  My  =  0.8.  Assume  ac  =  0. 


0.0015 


Equations  (3-19)  and  (3-20)  will  yield  a  value  for  the  drag  of  the  tail 
rotor.  Dy^  =  28.288. 

From  the  trim  values  obtained  earlier  and  using  equations  (3-21) 
through  (3-26)  determine  the  following  values  for  the  horizontal 
tai lplane. 


a  =  -0.0037  rad 

or y  =  -0.0349  rad 

C.  =  -0.1396 
LT 


cn  =  0.0216 

UT 

Ly  =  -12.7  lb 
Dy  =  19.6  lb 


Equations  (3-27)  through  (3-30)  will  yield  a  better  approximation  of  the 
main  rotor  lift  and  drag. 

Lp  =  10060  lbs  DF  =  -515.2  lbs 
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T 


Use  equations  (3-5)  through  (3-9)  to  obtain  better  extimates  for  the 
main  rotor  lift  and  drag  coefficients.  The  final  trim  values  of  the 
helicopter  are: 


M  =0.3  (-£-)=  0.053 


mt 

=  0.8 

(i, 

a 

=  -0.0026 

a 

=  0.092 

lf 

=  10060  lbs 

0! 

=  -8° 

df 

=  -515.2  lbs 

9.75 

II 

1 

e n 
o 

'V 

=  -0.045 

r  1  r  1 

The  trim  values  for  ,  — —  ,  and  9  are  used  to  enter  the  charts 
in  Section  7.5  of  Reference  5  to  get  the  nondimensional  isolated  deriv¬ 
atives  for  the  front  rotor. 


3p  J 


F 


-0.02 


3|J  J 


F 


0.014 


dk 


F 


-0.05 


0.17 


0.38 


C  1 

3<V 


3a 


c  J 


0.07 
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0.24 


Solidity  corrections  are  made  to  these  values,  if  necessary,  using 
equations  (3-31)  through  (3-41). 

Using  the  isolated  derivatives  just  obtained,  the  main  rotor  local 
derivatives  can  be  calculated.  Computation  is  done  by  equations  (3-42) 
through  (3-49).  The  following  results  will  be  obtained: 


-5.56 


lb-sec 

ft 


=  73666 


lb 

rad 


3.89 


lb-sec 

ft 


3Df  ik 

— ^  =  13570 
Oofp  rad 


3alF 

3^  =  °-00°244 


rad- sec 
ft 


0.24 


3M 


HUB 


=  17.54  lb-sec 


3MHUB 

=  17237 


ib-ft 

rad 


3a  1  p 

=  -0.09881 
3q 

Use  Figure  A-l  to  obtain  the  partial  derivative  values  of  C^,  C^,  and 
with  respect  to  or  by  obtaining  the  slope  of  the  line  at  the  trim  value 
of  ofpy<-  (-0.0312).  The  following  values  were  obtained: 
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ac 

da 


lFUS 

FUS 


0.0005  _  0. 0287 
degree  radian 


3C 


da 


-0.0005  =  -0.0287 
degree  radian 


ac 


aa 


mfus  _ 

FUS 


0.0058  _  0.3323 
degree  radian 


The  fuselage  local  derivatives  can  now  be  calculated  using  equations 
(3-50)  through  (3-55).  The  following  results  will  be  obtained: 


3L 

3u 


FUS  _ 
FUS 


-0.628 


lb-sec 

ft 


3Lfus 

daFUS 


280.8 


lb 

rad 


30 

3u 


FUS  _ 
FUS 


4.576 


lb-sec 

ft 


=  _84.23  lb 
aaFUS  rad 


8Mcus 

=  -36.21  lb-sec 

3ufus 


aM 


FUS 


3a 


FUS 


48848 


lb- ft 
rad 


Values  for  the  horizontal  tailplane  derivatives  are  calculated  using  the 
previously  obtained  values  for  lT  and  Dy  and  equations  (3-56)  through 
(3-59). 


=  -0.1251 


3LT  1h 

— I  =  11430  JiL 
3aT  J  rad 


aoT  ,l 

ao,  ft 


aoT  -lu 

=  -7.07 

3a,  raa 

1 
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The  tail  rotor  derivative  is  obtained  from  equation  (3-20a)  and  is 


30 


3u 


—  =  0.1545 


TR 


1 b-sec 
ft 


K 


FT 


=1.0  for  this  configuration  and  K^p 


=  0, 


therefore: 


3aFUS 

3a 


30(y 

~3a 


=  0.685 


3aFUS 

3u 


0.000513 


Equations  (3-64)  through  (3-85)  are  used  to  calculate  the  aircraft  total 
stability  derivatives.  Alternatively,  the  computer  programs  in  Appendix 
B  will  simplify  this  process.  XUXW  can  be  used  to  determine  the  X-force 
derivatives,  ZUZW  to  find  the  Z- force  derivatives,  and  MUWQ  to  find  the 
pitching  moment  derivatives.  Values  obtained  by  either  method  are  as 
fol lows: 


X  =  -8.6318  ■ b~|  - 

u  ft 


X  =  -19.075  1b'5ec- 

w  ft 


Z  =  4.344  4^ 

u  ft 


Z  =  -375.41  1b~^ec 
w  ft 


sec 


M  =  -4.783  ,. 

j  rt 


M  =  307.32 
w  ft 


M  -  -17838.2  -  C~ 
q  rad 
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MUWQ  incorporates  some  additional  relations  which  must  be  known  if 
equations  (3-80)  through  (3-85)  are  being  hand  calculated. 


8Mfus 

3u 


8MFUS  +  3MFUS 
3uFUS  acfFUS 


aaFUS 

3u 


3MHUB  3MHUB  9MHUB  3a 


3u 

3uF 

3aF 

3MFUS  . 

_  3Mfus 

aaFUS 

3a 

3o,fus 

3a 

3MHUB 

3Mhub 

3a 

F 

F 

F 

3a 

3a  F 

3a 

3MHUBf  _  3MHUBf  ^  x  3MHUBf  ebfi2Ms  ^da1^ 

3q  =  "  3U  1ZF  "  VT  ~da  \  +  T~  VStf, 


Divide  the  force  derivatives  by  mass  and  the  moment  derivatives  by 


pitching  moment  of  inertia  (I  )  to  form 

the  stability 

determi nant 

equation  (4-9). 

S  +  .0278 

0.0614 

32.2 

-0.014 

S  +  1.2079 

-203  S 

=  0 

0.0003 

-0.0176 

S2  +  1.019  S 

The  characteristic 

equation  will  be: 

S4  +  2.255 

S3  -  2.2788  S3  -  0.0776 

S  -  0.0037  =  0 
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Solving  the  quartic  equation  will  yield  the  following  roots: 


0.7843 

-3.0049 

-0.0172  ±  0.0357 


Therefore  this  helicopter  has  the  characteristics  of  a  heavily  damped 
root,  a  divergent  motion,  and  a  damped  oscillatory  motion. 

The  time  to  half  amplitude  of  the  damped  root  is: 

69 

\  ~  I -3700491  =  0,230  sec- 

The  time  to  double  amplitude  of  the  divergent  root  is: 


t2  = 


.69 

10.7843  f 


=  0.880  sec. 


The  time  to  half  amplitude  for  the  oscillatory  mode  is: 


.69 

1-0.01721 


=  40. 12  sec. 
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WORKSHEET  FOR  STAB 


WORKSHEET  FOR  XUXW 


APPENDIX  B 

COMPUTER  PROGRAMS  FOR  THE  HP-41CV  CALCULATOR 


These  computer  programs  are  written  for  the  HP-41  pocket  calculator 
and  are  intended  to  help  streamline  the  process  of  finding  solutions  to 
some  of  the  problems  generated  by  stability  and  control  analysis.  Input 
values  should  be  stored  as  indicated.  When  prompted  for  an  input  value, 
the  program  will  store  that  value  in  its  proper  storage  register.  When 
the  program  is  executed,  the  output  values  will  be  labeled  by  alpha 
characters. 

Worksheets  are  included  as  an  aid  to  organize  the  input  data  for  the 
programs  and  for  the  case  of  the  program  STAB  for  recording  output  data. 

A.  STAB  will  calculate  many  helicopter  parameters  needed  for  further 
calculations.  Both  main  rotor  and  tail  rotor  data  can  be  calculated. 
The  equations  solved  for  the  output  values  are  listed  in  the  definitions 
section  of  Chapter  III.  Input  values  and  their  storage  locations  are 
shown  in  the  worksheet  for  STAB.  Space  is  provided  on  this  worksheet  to 
record  the  output  values.  Output  values  are  labeled  by  alpha  characters. 

B.  XUXW  calculates  the  X-force  total  stability  derivatives.  Equations 
(3-64)  through  (3-71)  are  solved  with  this  program.  The  results  are 
displayed  using  alphanumerics  and  are  also  stored  for  use  in  other 
programs.  See  the  worksheet  for  XUXW  for  the  input  parameters  and  their 
proper  storage  location. 


C.  ZUZW  calculates  the  total  stability  derivatives  for  the  Z-direction. 
Input  data  is  the  same  as  that  for  XUXW  .  Output  data  is  listed  by 
alphanumerics  and  is  also  stored  for  future  use. 

D.  MUWQ  calculates  the  total  stability  derivatives  for  the  pitching 
moments  Mu,  and  M^.  A  worksheet  with  input  parameters  and  proper 
storage  locations  is  also  provided  for  this  program.  Output  data  is 
labeled  with  alphanumerics. 

E.  CE  finds  the  coefficients  of  the  characteristic  equation  generated 
by  the  stability  determinant,  equation  (4-9).  The  program  prompts  for 
the  input  values.  It  is  important  to  note  that  the  X  and  Z  derivatives 
must  be  normalized  by  aircraft  mass  ard  the  M  derivatives  must  be  nor¬ 
malized  by  Iy.  Also,  velocity  must  be  input  in  units  of  feet  per  second. 
The  output  data  will  be  the  coefficients  of  the  stability  quartic, 
equation  (4-9). 
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